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EXECUTIVE  SUMMARY 


The  purpose  of  this  project  was  to  assemble  a  test  gas  distribution  system  and  to 
perform  laboratory  experiments  with  an  existing  gas  phase  pulsed  streamer  corona  reactor  to 
support  the  US  Air  Force's  effort  to  investigate  the  feasibility  of  using  this  technology  for 
removing  nitrogen  oxides  from  the  combustion  exhaust  gases  of  jet  engine  testing  facilities. 
This  contract  was  performed  at  the  FAMU-FSU  College  of  Engineering  through  a  subcontract 
between  The  Florida  State  University  and  Applied  Research  Associates,  Inc.  (ARA).  This 
subcontract  was  created  from  a  primary  contract  between  ARA  and  the  United  States  Air 
Force.  This  project  is  a  continuation  of  Contract  FSESCR4 1580100,  Construction  of  a  Pulsed 
Streamer  Corona  Reactor,  performed  at  the  Department  of  Chemical  Engineering,  FAMU-FSU 
College  of  Engineering  (Locke  and  Finney,  1995). 

A  pulsed  streamer  corona  reactor  utilizes  a  high  voltage  electrical  discharge  in  a  non- 
uniform  electrode  geometry  to  initiate  chemical  reactions  that  lead  to  the  removal  of  various 
pollutants.  The  physical  aspects  of  a  gas  phase  pulsed  discharge  in  a  non-uniform  electrical 
field  include  the  formation  of  ionization  waves  (streamers)  through  the  growth  of  electron 
avalanches  formed  by  electron  impact  ionization  events  in  the  gas.  The  streamer  is  a  region  of 
highly  ionized  gas  (a  non-then-nal  plasma)  where  a  wide  range  of  highly  reactive  radicals  and 
species  are  formed  through  collisions  among  electrons,  molecules,  and  ions.  If  nitric  oxide 
(NO)  is  present  in  the  test  gas,  it  is  generally  first  converted  to  nitrogen  dioxide  (N02)  through 
a  variety  of  reactions  that  may  include  direct  reaction  with  dissociated  oxygen  or  reactions 
with  ozone.  N02may  be  removed  through  reactions  with  hydroxyl  radicals  created  from 
water  vapor  or  possibly  by  organic  species.  The  N02  is  converted  to  nitric  acid  aerosols  that 
in  turn  may  be  removed  by  scrubbers,  particle  filtration  devices,  or  electrostatic  precipitators. 

The  major  objectives  of  the  present  work  were:  1)  to  run  the  FSU-constructed  pulse 
streamer  corona  reactor  under  a  variety  of  gas  feed  conditions,  including  the  addition  of  water 
and  ethylene,  2)  to  determine  the  effect  of  the  presence  of  ethylene  on  nitrogen  oxide  removal. 
3)  to  look  for  byproduct  formation  from  the  reactions  of  ethylene  with  nitrogen  oxides,  and  4) 
to  develop  a  mathematical  model  of  the  gas  phase  pulsed  corona  reactor  that  accounts  for  the 
major  chemical  reactions  occurring  in  the  process.  In  addition,  basic  electrical 
characterization  of  the  corona  reactor  was  performed  in  order  to  determine  power  usage  and 
energy  efficiency  for  nitrogen  oxide  removal,  and  to  allow  comparison  with  other  data 
reported  in  the  literature. 

The  major  results  of  this  study  are  as  follows.  The  presence  of  500  ppm  ethylene  in 
dry  air  in  the  pulsed  corona  reactor  results  in  significant  enhancement  of  NO  removal. 
Quantitative  analysis  of  treated  gas  samples  by  GC/MS  indicates  that  23%  to  26%  of  the 
ethylene  is  decomposed  in  the  corona  reactor  with  an  input  voltage  of  50  kV  and  a  residence 
time  of  0.6  min.  The  presence  of  water  in  the  air  feed  (no  ethylene)  has  a  smaller  effect  than 
that  of  ethylene  in  dry  air,  and  the  presence  of  water  and  ethylene  in  air  shows  a  small 
improvement  in  NO  removal  than  the  case  with  ethylene  only  (no  water).  Ethylene 
breakdown  byproducts  were  not  found  by  sample  analysis  with  GC/MS.  Energy  costs  for  NO 


IV 


removal  for  dry  air,  air  with  ethylene,  air  with  water,  and  air  with  both  ethylene  and  water 
are  in  the  same  general  range  as  those  reported  in  the  literature.  The  mathematical  model 
developed  for  this  work  predicts  trends  similar  to  those  found  in  the  experiments  on  the 
effects  of  water  and  ethylene  on  NO  removal.  The  reaction  chemistry  in  the  model 
predicts  that  some  fraction  of  NO  is  converted  to  N2  and  02,  and  that  N02  levels  are 
reduced  in  the  presence  of  ethylene. 
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1.  GENERAL  INTRODUCTION  AND  BACKGROUND 


A  pulsed  streamer  corona  reactor  utilizes  a  high  voltage  electrical  discharge  produced 
within  a  non-uniform  electrode  geometry  to  initiate  chemical  reactions  that  lead  to  the 
removal  of  various  pollutants  from  a  gas  stream.  The  physical  aspects  of  a  gas  phase  pulsed 
discharge  in  a  non-uniform  electrical  field  include  the  formation  of  ionization  waves 
(streamers)  through  the  growth  of  electron  avalanches  formed  by  electron  impact  ionization 
events  in  the  gas  (Nassar,  1971;  Gallimberti,  1987;  Bastein  and  Marode,  1979,1995;  Morrow, 
1985).  The  streamer  is  a  region  of  highly  ionized  gas  (a  non-therrnal  plasma),  where  a  wide 
range  of  highly  reactive  radicals  and  chemical  species  are  formed  through  collisions  among 
electrons,  molecules,  and  ions.  NO  is  generally  first  converted  to  N02  through  a  variety  of 
reactions  that  may  include  direct  reaction  with  dissociated  oxygen  or  reactions  with  ozone. 
N02  may  be  removed  through  reactions  with  hydroxyl  radicals  created  from  water  vapor  or 
possibly  by  organic  species.  The  N02  is  thus  converted  to  nitric  acid  aerosols  that  in  turn  may 
be  removed  by  scrubbers,  particle  filtration  devices,  or  electrostatic  precipitators. 

The  physical  and  chemical  processes  involved  in  gas  phase  non-thermal  plasma 
processes  have  been  the  subject  of  a  number  of  experimental  and  theoretical  studies.  In  this 
section,  a  brief  review  of  the  literature  relevant  to  the  application  of  pulsed  streamer  corona  for 
the  removal  of  nitrogen  oxides  from  combusti  on  gases  will  be  presented.  The  basic  chemical 
reactions  expected  to  occur  in  the  process  and  a  description  of  the  physical  process  of 
streamer  propagation  will  be  given. 

1.1.  Chemical  Reactions  in  Pulsed  Corona 

The  chemical  reaction  mechanisms  for  the  removal  of  nitrogen  oxides  (NO  and  NO) 
by  gas  phase  non-thermal  plasma  reactions  have  been  studied  extensively  (Creyghton,  1994; 
Tas,  1995;  Alekseev  et  al.,  1993;  Wren,  1989;  McFarlane  and  Wren,  1991;  Mizuno  et  al., 
1993,  Tokanaga  and  Suzuki,  1984;  Matzing  1989).  One  model  of  the  reaction  chemistry  that 
has  been  extensively  used  for  pulsed  corona  systems  was  developed  at  IVTAN  in  Moscow 
(Alekseev  et  al.,  1993),  and  this  code  uses  over  900  individual  chemical  reactions.  The 
chemical  reactions  used  in  this  model  were  adapted  from  earlier  work  on  electron  beam 
technology  in  Germany  (Matzing,  1989).  The  extensive  nature  of  these  chemical  reactions 
arises  from  the  many  molecular,  ionic,  and  radical  species  that  are  formed  from  an  exhaust  gas 
that  may  include  S02,  NO,  N02,  CO,  C02,  02,  N2,  and  H20.  Recent  work  has  shown, 
however,  that  a  system  of  approximately  50  chemical  reactions  can  describe  the  overall 
reaction  very  well  (Creyghton.  1994;  Tas,  1995). 

Previous  modeling  and  experimental  work  with  pulsed  streamer  corona  reactors  and 
similar  non-thermal  plasma  systems  have  shown  that  nitrogen  oxide  removal  is  a  strong 
function  of  the  applied  electrical  field  and  the  composition  of  the  feed  gases.  The  amounts  of 
water  vapor  and  oxygen  in  the  feed,  as  well  as  the  presence  of  some  organic  species,  have 
been  found  to  be  very  important  in  affecting  the  nitrogen  oxide  removal  rates.  An  essential 
feature  of  any  method  that  will  improve  the  operation  of  the  pulsed  streamer  corona  reactor 
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for  nitric  oxide  removal  is  the  suppression  of  any  back  reactions  or  any  radical  recombination 
reactions  that  lead  to  the  formation  of  NO. 


High  energy  electrons  produced  in  the  pulsed  discharge  will  collide  with  the  major  gas 
phase  species,  primarily  N2  and  02  to  lead  to  dissociation  of  these  species  by 


Ni  +  e — >  2N  +  e- 

0) 

02  +  e-  20  +  e- 

(2) 

These  dissociated  nitrogen  and  oxygen  radicals  can  recombine,  react  together,  or  react  with 
other  species. 

The  primary  reactions  that  lead  to  nitrogen  oxide  removal  consist  first  of  the  oxidation 


of  NO  to  N02  through  reactions  with  02,  03,  and  O. 

2NO  +  02  -»  2N02  (3) 

NO +0  +  M  ->  N02 +M  (4) 

NO  +  03  ->•  NOj  +  02  (5) 

In  the  presence  of  water  vapor,  NO  can  also  be  oxidized  by  the  hydroperoxy  radical  through 

NO  +  HO,'  ->  N02  +  OH'  (6) 


and  subsequent  reactions  of  N02  with  hydroxyl  radical  lead  to  acid  formation 

N02  +  OH'  -»  HNOj  (7) 

The  nitric  acid  thus  formed  may  condense  in  existing  water  droplets,  it  may  form  acid 
aerosols  through  heterogeneous  nucleation  via  existing  particles  in  the  gas  (Kiang  et  al., 
1973),  or  it  may  form  acid  aerosols  via  homogeneous  nucleation  and  particle  growth 
(Seinfeld,  1986). 

In  the  presence  of  water  vapor  and  oxygen,  direct  electron  impact  with  molecular 

oxygen 

e-  +02  »  0  +  0(’D),  (8) 

and  subsequent  reaction  of  singlet  oxygen  with  water  vapor  leads  to  hydroxyl  radical 
production 

0('D)  +  H20  ->  20H'.  (9) 
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Hydroperoxy  radicals  can  also  be  formed  though  reactions  of  hydroxyl  radicals  and  carbon 
monoxide 


CO  +  OH'  +  02  ->  C02  +  H02. 


(10) 


The  oxidation  of  most  hydrocarbons  also  lead  to  formaldehyde  (Seinfeld,  1986)  which  can  in 
turn  lead  to  additional  hydroperoxy  radicals  via 


hcho  +  oh'  ho;  +  co  +  e2o.  (ii) 

Many  other  reactions  may  be  occurring  in  the  system  depending  upon  the  composition  of  the 
feed  gas  (Matzing,  1989;  Alekseev  et  al.,  1993;  McFarlane  and  Wren,  1991),  and  of  particular 
importance  are  the  recombination  and  back  reactions  that  produce  additional  NO  that  will 
tend  to  reduce  the  efficiency  of  the  process. 

In  dielectric  barrier  discharge  reactors  (Chang  et  al.,  1992),  and  in  pulsed  corona 
reactors  (Masuda  and  Nakao,  1986;  Shimazu  and  Mizuno,  1993),  it  has  been  established  that 
the  addition  of  water  vapor  will  greatly  improve  the  removal  rates  of  NO.  This  is  considered 
to  be  primarily  due  to  the  enhanced  production  of  hydroxyl  radicals.  Optimal  values  of  water 
vapor  will  depend  upon  both  temperature  and  electric  field  conditions.  Increasing  02  content 
has  also  shown  improvement  in  NO  removal.  In  dielectric  baffler  discharges,  increases  in  02 
content  show  an  optimal  removal  near  5%  (volume/volume)  of  oxygen  at  higher  applied 
voltages  for  dry  gas.  This  is  thought  to  be  due  to  the  increased  production  of  hydroxyl 
radicals  through  reaction  with  singlet  oxygen  as  shown  above  (9);  however,  at  very  high 
oxygen  levels  generation  of  NO  will  occur  through  back  reactions  such  as 

N02  +  0  ->  NO  +  02,  (12) 

or  recombination  reactions  such  as 


N  +  0  NO, 

(13) 

both  of  which  will  lead  to  a  reduction  in  efficiency  of  NO  removal.  The  above  two  reactions 
show  that  it  is  desirable  to  remove  N02  as  soon  as  it  is  formed,  and  that  optimal  levels  of  O 
must  be  achieved  to  allow  hydroxyl  radicals  to  be  formed  through  reactions  of  0  with  water, 
but  to  not  allow  too  much  NO  to  be  formed.  NO  can  also  be  formed  through  such  reactions  as 

N2  +  0  -*  NO  +  0 

(14) 

N  +  02  -*  NO  +  0 

05) 

N  +  OB'  -»  NO  ■+•  H, 

(16) 
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which  is  the  well  known  Zeldovich  mechanism  (Zeldovich,  1946)  developed  for  thermal 
combustion  devices,  and  therefore  dissociated  N  as  well  as  O  can  lead  to  additional  NO 
formation. 

It  is  also  interesting  to  note  that  N02  can  react  with  N  (McFarlane  and  Wren,  1991)  to 
regenerate  N2  by 

N02  +  N  -+  N2  +  20  (17) 

N02  +  N  -*■  N2  +  02.  (18) 

Due  to  the  somewhat  slower  rate  constants  for  these  reactions,  they  are  generally  more 
important  at  longer  residence  times.  The  enhancement  of  these  reactions  is  certainly  desirable 
because  of  the  obvious  advantages  of  the  final  products.  Tas  (1995)  concluded  that  reduction 
of  NO  did  not  occur  in  his  system  because  he  could  not  detect  N2  production.  He  performed 
experiments  in  He  and  N2  environments,  however,  it  is  unlikely  that  N2  would  be  produced  in 
the  He  environment,  and  that  in  the  N2  environment  the  addition  of  even  several  hundred  parts 
per  million  N2  would  be  extremely  difficult  to  detect.  The  question  of  how  much,  if  any,  N2  is 
produced  by  the  above  reactions  remains  open. 

The  oxidation  of  NO  to  N02  has  been  found  to  occur  only  after  the  onset  of  the  corona 
discharge  (Masuda  and  Nakao,  1990).  Increases  in  peak  voltages  beyond  corona  onset  were 
found  to  maintain  total  NOx  at  a  constant  until  a  critical  field  intensity  was  reached.  This 
indicates  that  NO  was  oxidized  to  N02  without  further  conversion  of  N02  to  acid.  Beyond 
this  critical  field  intensity  the  N02  value  rapidly  decreased  (Masuda  and  Nakao,  1990).  These 
authors  claim  that  the  highest  possible  peak  voltage  will  result  in  the  most  effective  removal 
strategy  in  order  to  minimize  the  amount  of  NO  formed  by  recombination  of  N02  and  O; 
however  the  use  of  water  and  or  ammonia  is  necessary  to  scavenge  the  N02  before  it  can  react 
with  oxygen  to  regenerate  NO.  One  would  expect  that  at  very  high  fields  NO  will  be  formed 
through  the  recombination  of  N  and  O  and  that  there  must  therefore  be  an  optimal  field 
strength.  Indeed,  Ohkubo  et  al.  (1994)  observed  an  optimum  in  NOx  reduction  as  a  function 
of  corona  power  for  certain  treatment  conditions,  and  they  also  observed  a  general  decrease  in 
energy  yield  of  NO  removal  with  corona  power.  By  optimizing  the  addition  of  NH3,  C02,  02, 
and  N2  they  were  able  to  show  increased  NO  removal  and  higher  energy  yields. 

Recent  work  by  the  group  at  Eindhoven  Technical  University  (Tas,  1995;  Creyghton. 
1994)  have  considered  the  effects  of  water  vapor,  oxygen  content  and  the  addition  of  solid 
catalytic  particles  on  NO  removal  in  pulsed  corona  reactors.  They  found  that  increasing  the 
oxygen  content  from  0%  to  9%  caused  a  decrease  in  efficiency  of  NO  removal  in  a  dry 
nitrogen  feed  gas.  When  water  vapor  was  added,  they  found  the  opposite  trend;  the  removal 
efficiency  increased  as  the  water  content  was  varied  from  0.0  to  11.7%.  They  did  not  find 
optimal  values  for  02  or  H20  as  was  observed  in  the  work  on  dielectric  discharge  reactors. 
Their  corona  reactor  was  run  with  constant  residence  time  (1  see)  and  pulsed  voltage  level  (28 
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kV),  constant  NO  input  (1200  ppm),  and  with  variable  pulsed  voltage  frequency  (5  -  240  Hz). 
The  power  input  to  the  reactor  varied  from  0.05  to  1.5  watts. 

A  number  of  studies  have  emphasized  the  importance  of  the  presence  of  ammonia  to 
enhance  the  removal  of  nitrogen  oxides  via  reactions  with  nitric  acid  to  form  ammonium 
nitrate 

HNOj  +  NHj  ->  NH4NO,  (19) 


which  will  create  a  solid  aerosol  particle  that  can  be  removed  by  electrostatic  precipitation  or 
filtration.  The  advantage  of  this  process  is  the  production  of  a  salable  byproduct  that  can  be 
used  as  a  fertilizer;  however  the  disadvantages  are:  1)  the  need  for  handling  large  quantities  of 
highly  toxic  ammonia  gas,  and  2)  the  necessity  to  remove  any  unreacted  ammonia  from  the 
exit  gas  stream.  Since  the  addition  of  ammonia  is  not  desired  in  the  present  work  due  to  these 
disadvantages,  further  analysis  of  the  addition  of  water  vapor  and  oxygen  will  be  considered. 
It  is  also  important  to  note  that  the  work  by  Masuda  and  Nakao  (1990)  did  not  consider  the 
role  of  CO  which,  as  shown  above,  is  very  important  to  the  formation  of  hydroperoxy 
radicals. 

It  is  clear  that,  depending  upon  the  composition  of  the  feed  and  the  reactor  conditions, 
there  are  optimal  levels  of  such  additives  as  water  vapor,  oxygen,  and  possibly  hydrocarbons 
(or  CO)  that  will  most  effectively  promote  the  removal  of  nitrogen  oxides  as  acid  aerosols.  It 
is  also  clear  that  the  role  of  hydrocarbons,  through  the  formation  of  formaldehyde  and  carbon 
monoxide,  will  enhance  the  production  of  hydroperoxy  radicals  through  reactions  with 
hydroxyl  radicals.  In  addition,  the  electric  field  conditions  in  the  reactor  will  also  strongly 
affect  the  removal  of  the  nitrogen  oxides  as  shown  in  the  work  of  Masuda  and  Nakao  (1990) 
and  Mukkavilli  et  al.  (1988).  It  is  also  clear  that  in  order  to  enhance  the  efficiency  of  the 
removal  of  NO  it  is  necessary  to  suppress  the  "back"  reactions  and  thus  minimize  the 
production  of  additional  undesired  NO.  This  can  be  accomplished  through  reactions  of  N02 
with  hydroxyl  radicals,  and  the  use  of  optimal  levels  of  oxygen.  Precise  control  of  water 
vapor  and  oxygen  input  would  thus  serve  as  a  means  to  optimize  the  removal  of  NO  and 
enhancement  of  reactions  that  lead  to  N  and  O  products  is  also  desirable. 

The  chemistry  of  ethylene  has  been  studied  in  relation  to  gas  phase  atmospheric 
chemistry  (Atkenson  and  Lloyd,  1984;  Seinfeld,  1986).  The  primary  reactions  of  ethylene 
with  hydroxyl  radicals,  ozone,  and  dissociated  oxygen  and  the  resulting  effects  on 
photochemical  smog  have  been  determined.  To  date  only  one  report  by  Mizuno  (Mizuno  et 
al.,  1993)  has  considered  the  effects  of  a  pulsed  corona  discharge  on  hydrocarbon  reactions. 
Some  work  has  been  reported  on  VOC  removal  with  pulsed  corona  technology  (Futamura  et 
al.,  1995). 

In  the  natural  atmosphere,  ethylene  is  expected  to  react  primarily  with  hydroxyl 
radicals  since  the  concentration  of  singlet  oxygen  and  ozone  are  generally  small.  The  reaction 
pathway  for  ethylene  in  the  presence  of  oxygen,  NO  and  hydroxyl  radicals  has  been  well 
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studied  (Seinfeld,  1986).  Hydroxyl  radicals  directly  attack  ethylene  to  produce  peroxyl 
radicals  that  can  further  react  with  oxygen  and  NO  through  the  following  reactions 


C2H4  +  OH'  ->  HOCH2CH2  (20) 

HOCH,CH2'  +  02  ->  H0CH2CH202.  (2l) 

NO  is  converted  to  N02  through 

hoch2ch2o2’  +  no  ->  no2  +  hoch2ch2o:  (22) 

Approximately  72%  of  the  ethylene  leads  to  formaldehyde  production 

H0CH2CH20'  +  02  ->  HCHO  +  CH2OH;  (23) 

while  the  remaining  28%  leads  to 

H0CH2CH20'  +  02  HOCH2CHO  +  HOj  (24) 

CH2OH'  +  02  HCHO  4-  H02.  (25) 

The  formaldehyde  produced  in  the  above  reactions  can  further  react  with  hydroxyl  radicals  to 
produce  hydroperoxy  radicals  as  discussed  previously  (see  Equation  11). 

Direct  reactions  of  ethylene  with  ozone  are  also  known  to  occur.  These  reactions  lead 
to  the  production  of  formaldehyde,  formic  acid,  carbon  dioxide,  and  carbon  monoxide  by 

C2H4  +  Os  ->  HCHO  +  HXOO  (26) 

H2COO  ->  CO*  +  H2  (27) 

H2COO  ->  CO  +  H20  (28) 

H2COO  +  20,  CO,  +  2HO,,  (29) 

with  a  net  reaction  given  by 

C2H4  +  O,  ->  HCHO  +  0.4H2COO  +  0.18CO2 

+  0.42CO  4-  1.2H2  +  0.42H,O  4-  0.12HO2.  (30) 

If  no  water  is  initially  present  in  the  feed,  the  water  formed  through  the  ozone  reactions  with 
ethylene  may  lead  to  the  production  of  hydroxyl  radicals. 
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Mizuno  et  al.  (1993),  performed  pulsed  corona  reactor  experiments  with  dry  air 
containing  NO  and  ethylene.  He  observed  the  formation  of  acetic  acid,  however,  none  of  the 
other  products  considered  above  were  directly  observed  using  Fourier  IR  Spectroscopy. 
Additional  reactions,  normally  considered  negligible  in  the  atmosphere,  include 


C2H4  +  0(3P)  ->  CH3  +  HCO 

(31) 

C2H4  +  0(3P)  -►  CH2CHO  fH 

(32) 

C2H4  +  0(3P)  CH2CO  +  H, 

(33) 

These  reactions  may  be  important  in  the  plasma  environment  of  the  pulsed  corona  reactor. 
According  to  Atkenson  and  Lloyd  (1984),  the  reaction  pathways  of  the  above  reactions  gives 
55%  to  the  first,  36%  to  the  second,  and  the  remainder  to  the  third  reaction.  The  vinoxy 
radical  at  low  pressure  leads  to 

CH2CHG  +  02  ECHO  +  CO  +  OH  (34) 

1.2.  Physical  Processes  Occurring  in  Streamer  Propagation 

The  chemical  reactions  initiated  by  a  pulsed  corona  discharge  that  were  discussed 
above  are  a  direct  result  of  a  number  of  physical  processes  that  occur  in  the  non-thermal 
plasma.  The  purpose  of  this  section  is  to  provide  A  brief  review  of  some  of  the  basic  physical 
processes  occurring  in  a  pulsed  corona  reactor.  The  material  is  obtained  primarily  from 
Nasser  (1971)  and  Creyghton  (1994). 

Normal  air  at  standard  conditions  contains  approximately  103  electrons,  positive  ions, 
and  negative  ions  per  cm3.  These  species  arise  from  UV  and  cosmic  radiation  and  are 
generally  in  a  Boltzmann  equilibrium.  In  the  presence  of  a  small  uniform  DC  electric  field  (1 
V/cm),  small  current  pulses  of  less  than  10'16  A  arise.  When  enough  free  electrons  are  present 
due  to  external  volume  ionization  and  electrode  emissions  from  external  radiation,  a  steady 
current  arises.  In  this  situation,  the  current  flow  due  to  the  motion  of  charged  species  has  a 
small  effect  on  the  Boltzmann  equilibrium,  and  the  conductivity  depends  upon  the  rate  of  ion 
and  electron  production,  recombination,  and  motion.  As  the  applied  electrical  field  increases, 
the  equilibrium  is  upset  as  the  number  of  ions  and  electrons  are  neutralized  when  they  reach 
the  electrodes.  The  rate  of  increase  of  current  with  voltage  decreases  and  a  saturation  region 
with  current  density  of  approximately  10'9  A/cm2  occurs.  This  is  termed  the  "dark  discharge 
region",  which  is  a  region  where  the  discharge  is  not  self-sustaining  because  it  requires  an 
external  source  of  radiation  to  provide  electrons  and  ions. 

A  further  increase  in  voltage  causes  the  current  to  rise  at  an  increasing  rate.  Ionization 
by  collision  or  electron  impact  ionization  occurs  through 

A  +  e-  A4  +  c*  +  e\  (35) 
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whereby  an  avalanche  is  created  through  a  cascade  where  two  electrons  are  released  for  every 
one  impact  with  molecular  species  A.  This  is  known  as  the  Townsend  discharge  region  and 
the  current  is  generally  greater  than  1  pA.  This  region  is  followed  by  an  abrupt  transition 
leading  to  breakdown,  where  an  over  exponential  increase  in  current  is  followed  by  a  collapse 
in  voltage.  The  current  can  be  increased  by  reducing  the  resistors  of  the  outer  circuit  which 
will  lead  to  a  voltage  drop  and  a  glow  region.  Higher  currents  will  lead  to  arc  formation. 

A  number  of  different  types  of  corona  discharges  can  occur  in  non-uniform  electric 
fields,  and  these  are  reviewed  in  the  literature  (Nasser,  1971).  Of  direct  interest  here  is  pulsed 
corona,  which  is  characterized  by  strong  branching  and  a  large  region  of  existence  in  terms  of 
the  applied  voltage. 

The  key  aspect  to  understanding  a  corona  discharge  is  a  proper  description  of  the 
Townsend  avalanche.  Yamamoto  states  "...the  Townsend  avalanche  is  a  conceptual 
description  of  the  electron  multiplication  process.  The  process  begins  with  a  free  electron  in  a 
strong  field.  The  electron  gains  enough  energy  through  collisions  with  gas  molecules  to 
ionize  one  of  them.  This  produces  another  electron  and  a  positive  ion.  The  two  free  electrons 

are  able  to  produce  further  electrons  in  a  cascade . The  rate  of  ionization  is  described  by  the 

first  Townsend  coefficient  (X,  which  is  the  number  of  electrons  produced  per  cm'  per  initial 
electron."  The  major  factors  that  limit  the  growth  of  avalanches  include  1)  a  rapid  decrease  in 
field  as  the  distance  from  the  electrode  increases,  2)  electron  attachment  to  gas  molecules 
which  removes  electrons  from  the  avalanche  leading  to  the  repulsion  of  electrons  from 
negative  ions,  and  3)  a  buildup  of  positive  ions  behind  the  front  of  electrons  which  will  lead  to 
a  retardation  of  the  rate  of  electron  propagation  (Nasser,  1971). 

The  transition  from  an  avalanche  to  a  "streamer"  in  a  uniform  electric  field  begins  at 
the  cathode  where  the  avalanche  is  initiated.  This  avalanche  propagates  to  the  anode,  and  near 
the  anode  a  build-up  of  positive  ions  occurs.  During  the  primary  avalanche,  propagation 
excitation  of  gas  molecules  occurs.  Photons  are  emitted  by  these  excited  states  and  these 
photons  lead  to  photoionization  and  the  creation  of  photo  electrons.  These  in  turn  create 
secondary  avalanches  that  will  also  propagate  toward  the  anode  and  lead  to  additional 
photo  ionization  and  third  generation  avalanches.  The  merging  of  the  secondary  avalanches 
with  positive  space  charge  gives  rise  to  the  growth  of  an  ionized  channel  from  the  anode  to  the 
cathode,  and  the  electrons  w  ill  move  toward  the  anode.  The  propagation  of  one  streamer  tip 
continues  and  the  others  stop  advancing  due  to  the  lack  of  avalanches  feeding  them. 

Of  the  mechanisms  for  ionization,  photoionization,  given  by 

A  +  hv  -»  A  +  e  (36) 

is  the  most  important.  The  electric  field  generated  by  the  negative  ions  allows  streamers  to 
propagate  into  regions  of  low  field.  The  electron  avalanche  forms  streamers  when  the  electric 
field  created  by  the  electron  space  charge  reaches  the  magnitude  of  the  geometric  field.  In  the 
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words  of  Nasser  (1971),  "...the  streamer  discharges  are  highly  localized  positive  or  negative 
space  charge  waves  which  enhance  the  applied  field  in  front  of  the  wave  (active  region)  and 

propagate  because  of  electron  avalanches  in  the  high  field . photoionization  is  generally 

assumed  to  be  the  mechanism  that  supplies  secondary  seed  electrons  and  triggers  avalanches. 
Along  the  track  of  the  wave  remains  a  weakly  ionized  filament  plasma  (passive  region), 
along  which  the  conduction  current  flows  to  the  high  voltage  electrode,  and  supplies 
energy  for  streamer  advancement." 

The  major  characteristics  of  streamers  include  length,  number,  speed  of 
propagation,  width,  number  of  electrons  and  ions,  time  lag,  and  electric  field  distribution. 
Typical  values  of  streamer  width  are  in  the  range  of  25  to  100  Em  and  the  speed  of 
propagation  is  in  the  range  of  0.7  -  0.9  mm/ns  (Creyghton,  1994).  The  number  of 
electrons  and  ions  varies  in  space  and  time  in  the  streamer  with  an  upper  range  of  about 
102O/m3  .  The  distribution  of  streamers  along  a  wire  electrode  has  been  measured  using 
Schlieren  photography  to  be  approximately  7/cm  (Creyghton,  1994).  The  Raether-Meek 
criterion  for  streamer  onset  in  a  uniform  field  has  been  extended  for  non-uniform  fields 
that  are  required  to  produce  pulsed  streamer  corona  (Creyghton,  1994).  This  onset  voltage 
for  a  non-uniform  field  in  a  wire-cylinder  geometry  according  to  this  criterion  is  about  17 
kV  (Creyghton,  1994). 

Creyghton  (1994)  reviews  the  basic  physical  models  of  streamer  formation.  In 
general,  one  dimensional  convective/diffusive/reaction  equations  for  positive  and  negative 
ions  and  electrons,  coupled  with  the  Poisson  equation  for  the  electric  field  and  given 
electron  velocity  distributions  (obtained  experimentally  by  solution  of  the  Boltzmann 
equation,  or  by  Monte  Carlo  simulation)  are  solved  to  determine  the  rate  of  propagation  of 
the  streamer  in  a  uniform  or  non-uniform  electric  field.  Simulations  have  been  performed 
coupling  streamer  models  with  chemical  kinetics  (Alekseev  et  al.,  1993,  McFarlane  and 
Wren,  1991).  Due  to  a  number  of  uncertainties  in  the  properties  of  the  streamers,  no 
direct  model/experimental  comparisons  have  been  made  without  some  decree  of  empirical 
fitting. 
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2.  EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


2.1.  Power  Supply  and  Corona  Reactor 

The  reactor  and  power  supply  used  in  the  present  study  was  constructed  on  contract 
with  the  United  States  Air  Force  and  is  described  in  extensive  detail  by  Locke  and  Finney 
(1995).  The  overall  experimental  system  is  shown  in  Figure  1.  Feed  gas  is  prepared  using  a 
building  supply  of  dry  air  and  gas  regulators/cylinders.  The  gases  are  metered  using  mass 
flow  controllers  and  then  are  fed  into  a  mixing  chamber.  Thereafter  the  feed  gas  flows  into 
the  corona  reactor.  Downstream  from  the  reactor  some  gas  is  vented  to  a  fume  hood,  some 
gas  flows  into  a  NOx  monitor,  and  some  gas  is  collected  in  sampling  vessels. 

The  pulsed  power  supply,  shown  schematically  in  Figures  2  and  3,  utilizes  a  high 
voltage  AC  60  Hz  input  provided  by  a  Universal  Voltronics  130  kV,  28  mA  power  supply. 
The  AC  voltage  first  flows  through  a  bank  of  333  kohm  current  limiting  resistors.  A  high 
voltage  rectifier,  made  up  of  a  diode  chain,  serves  to  modify  the  AC  voltage  waveform  by 
removing  the  lower  half  of  the  sine  wave.  A  capacitor  bank  located  between  line  voltage  and 
ground  charges  during  the  rising  portion  of  each  voltage  waveform.  A  mechanical  rotating 
spark  gap  matched  in  frequency  to  the  line  voltage  (60  Hz)  serves  as  the  pulse  producing 
device.  When  the  rotating  arm  and  the  sphere  electrodes  are  not  aligned,  the  voltage  charges 
the  capacitors;  the  gap  discharges  the  capacitors  when  the  arm  and  spheres  are  aligned.  The 
characteristics  of  the  capacitor  bank  determine  the  rise  time  for  the  pulsed  waveform  that  is 
ultimately  delivered  to  the  reactor.  The  rise  time  is  approximately  20  -  50  ns,  and  the  pulse 
duration  is  on  the  order  of  200  ns.  The  spark  gap  is  aligned  and  synchronized  to  line  voltage 
with  a  strobe  lamp.  The  unit  developed  in  this  work  is  fundamentally  similar  to  other  units 
developed  at  Florida  State  University  (Clements  et  al.,  1989;  Mizuno  and  Clements,  1987)  as 
well  as  to  those  used  in  the  literature  (Creyghton,  1994;  Tas,  1995;  Mizuno  et  al.,  1993). 

The  main  body  of  the  pulsed  corona  reactor  is  shown  in  Figure  4.  The  reactor  is 
constructed  of  a  4"  diameter,  #316  stainless  steel  cylinder  18"  long  from  end  to  end.  The 
active  region  of  the  reactor  consists  of  a  12"  section  of  the  reactor  where  a  stainless  steel  wire 
electrode  runs  concentrically  down  the  center  of  the  reactor  body.  Three  flow  rates  were  used 
in  the  present  study:  500  seem  (standard  cubic  centimeters  per  minute),  5000  seem,  and 
10,000  seem.  The  corresponding  residence  times  in  the  active  region  of  the  reactor  are  12 
min,.  1.2  min,  and  0.6  min  (36  seconds).  Corresponding  Reynolds  numbers  in  the  reactor  are 
7,  70,  and  140,  respectively.  This  indicates  that  for  all  runs  the  reactor  flow  is  laminar. 

2.2.  Gas  Feed  System 

The  gas  feed  system  consists  of:  1)  building  supply  source  of  dry  compressed  air,  2) 
gas  cylinders  and  regulators  containing  1%  NO  in  dry  N2  and  pure  ethylene,  3)  gas  cylinders 
and  regulators  containing  a  calibration  gas  of  100  ppm  NO  in  N2,  4)  water  tank  through  which 
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the  air  and  other  gas  can  be  bubbled,  5)  flow  meters  and  mass  flow  controllers,  and  6) 
mixing  chamber.  MKS,  Instruments,  Inc.,  mass  flow  controllers  are  used  on  the  bulk  air 
flow,  the  NO  gas  feed,  and  the  ethylene  gas  feed.  Bypass  loops  around  the  MKS  units  are 
available  if  needed.  The  mixing  chamber  contains  stainless  steel  wool  to  provide  a 
turbulent  environment  for  adequate  gas  mixing.  Figure  5  shows  a  schematic  of  the 
complete  gas  feed  system. 

2.3.  Analytical  Instrumentation 

The  analysis  of  the  outlet  gases  from  the  reactor  (Figure  5)  include  direct  measurement 
of  NO,  03,  N20,  and  hydrocarbon  breakdown  products.  Nitrogen  oxides  are  measured  using  a 
Thermo-Environmental  Instruments,  Inc.  Model  42  H  Chemiluminescence  NO  Analyzer. 
Ozone  is  measured  using  a  PCI  Ozone  Corporation,  Model  HC-1  ozone  monitor.  The  gas 
phase  products,  including  N20  and  hydrocarbon  breakdown  products,  are  measured  using  a 
Hewlett-Packard  GC/MS  instrument  located  in  the  Department  of  Chemistry  at  Florida  State 
University.  The  pulse  waveform  characteristics  are  measured  using  a  Tektronix  TDX  460 
Four  Channel  Digitizing  Oscilloscope  with  a  P6015A  lOOOx  3.0  pF,  1000  MO  voltage  probe 
and  a  TM502A/AM503B  current  probe.  Graphic  devices  included  a  Tektronix  HC-100  color 
plotter  for  printing  oscilloscope  waveforms,  and  a  Fisher  Recordall  Series  5000  chart  recorder 
to  monitor  NO  concentrations  from  the  NOx  monitor. 

2.4.  Sampling  Methods 

The  outlet  gas  from  the  reactor  was  fed  directly  into  the  Thermo-Environmental 
Instruments  NOx  analyzer  where  NO  and  N02  were  measured.  Gas  samples  were 
collected  in  500  ml  glass  sampling  chambers  equipped  with  ports  for  GC/MS  analysis. 
These  samples  were  taken  directly  to  the  Chemistry  Department  for  analysis  and  were 
usually  measured  within  several  hours  of  collection.  Controls  were  run  to  ensure  that 
leakage  and  degradation  within  the  sampling  containers  were  minimal.  The  sampling 
containers  were  washed  with  deionized  water  and  dried  overnight  with  clean  air. 

2.5.  Experimental  Run  Procedures 

1.  Turn  on  mass  flow  controllers  and  let  warm  up  for  at  least  30  minutes. 

2.  Open  the  valve  for  the  compressed  air  and  adjust  the  pressure  regulator  to  about 
20  psig. 

3.  Adjust  needle  valves  (2)  in  the  compressed  air  line  to  make  sure  air  at 
atmospheric  pressure  is  going  through  the  NOx  monitor. 

4.  Turn  on  the  NOx  monitor  and  let  warm  up  for  around  2  hours. 

5.  Attach  the  gas  calibration  standards  (NO,  N02)  to  the  NOx  monitor. 
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6.  Calibrate  the  N0X  monitor  to  the  concentrations  of  the  calibration  gas  (100 
ppm). 

7.  Turn  on  the  chart  recorder  and  adjust  to  necessary  paper  speed. 

8.  Once  the  mass  flow  controllers  stabilize,  make  sure  they  are  reading  zero  and 
adjust  as  necessary. 

9.  Clean  the  electrodes  in  the  rotating  spark  gap  using  sandpaper  and  ethanol. 

10.  Turn  on  the  rotating  spark  gap,  the  strobe  lamp,  the  fan  inside  the  faraday  cage, 
the  fan  on  the  diode  board,  and  the  fan  on  the  load  resistors. 

11.  Adjust  the  rotating  spark  gap  and  align  in  a  position  perpendicular  to  the 
discharge  position  when  strobe  lamp  is  in  line  frequency  mode  (60  Hz). 

12.  Check  that  the  valves  in  the  sample  gas  line  are  open,  and  make  sure  that  the  gas 
has  an  open  path  through  either  the  rotameters  or  the  mass  flow  controllers, 
through  the  mixing  chamber  and  the  reactor,  and  then  to  a  vent  or  into  an 
analytical  or  sampling  device. 

13.  On  the  gas  cylinder  that  is  to  be  used,  open  the  main  valve  on  the  tank. 

14.  Adjust  the  pressure  regulators  and  needle  valves  to  desired  flow  rate. 

15.  Open  sampling  vessel  and  let  sample  gas  pass  through  for  enough  time  to  come 
to  equilibrium. 

16.  Close  the  valves  on  the  sampling  tube  and  replace  with  another  sampling  tube. 

1 7.  Get  the  baseline  reading  from  NOx  monitor. 

18.  Turn  on  the  oscilloscope,  place  the  voltage  probe  at  the  insulator  conductor  on 
the  reactor  opposite  the  input  of  the  pulsed  power,  and  place  the  current  probe 
on  the  wire  grounding  the  reactor. 

19.  Ground  the  reactor  body  and  check  all  of  the  gas  and  electrical  connections 
between  the  pulse  generator  and  the  reactor. 

20.  Plug  in  the  AC  line  cord  for  the  power  supply. 
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21.  Put  all  three  switches  on  the  front  of  the  power  supply  in  the  ON  position, 
press  the  HV  reset  button,  press  the  HV  ON  button,  and  then  turn  the  dial  to  the 
desired  voltage. 

22.  Allow  all  analytical  devices  to  reach  equilibrium  for  each  new  condition  during 
experimental  run. 

23.  Adjust  waveforms  on  oscilloscope  and  save  the  waveforms  in  memory. 

24.  Close  the  valves  on  the  sampling  vessel  and  remove  from  the  system. 

25.  Turn  off  the  power  to  the  reactor. 

26.  Ground  the  reactor  and  rotating  spark  gap  with  the  portable  grounding  rod. 

27.  Plug  in  printer  and  attach  to  oscilloscope. 

28.  Print  hardcopy  of  waveforms. 

29.  Close  the  vanes  on  the  gas  cylinder. 

30.  Purge  the  system  with  compressed  air  for  several  minutes. 

3 1 .  Turn  off  the  mass  flow  controllers  and  the  NOx  monitor. 

32.  Close  the  compressed  air  valve  and  reduce  regulator  pressure. 
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3.  EXPERIMENTAL  RESULTS 


3.1.  Initial  Power  Supply  Characterization 

Preliminary  characterization  of  the  pulsed  voltage  waveform  supplied  to  the  corona 
reactor  was  carried  out  by  Cadet  Steve  Lucas  from  the  Department  of  Chemistry,  United 
States  Air  Force  Academy,  while  he  was  working  at  Florida  State  University  as  part  of  the 
Cadet  Summer  Research  Program  during  the  time  of  this  project.  Additional  work  including 
measurement  of  the  power  input  to  the  reactor  has  been  performed  recently,  and  all  of  the 
major  features  of  the  pulsed  power  supply  characterization  are  reported  below. 

The  initial  characterization  of  the  pulsed  power  supply  included  calibration  of  the 
power  supply  dial  settings,  and  measurement  of  the  voltage  pulse  width  (decay  time),  pulse 
rise  time,  resonant  frequency,  and  output  voltage  at  several  locations  in  the  circuit  and  for  a 
range  of  input  voltage  levels. 

Calibration  of  the  power  supply  dial  settings  showed  that  the  average  peak  AC  voltage 
delivered  to  the  pulse-forming  network  was  linear  with  the  dial  settings.  This  is  shown  in 
Figure  6  over  the  range  of  5  kV  to  30  kV.  Above  30  kV  the  average  peak  voltage  was 
determined  through  extrapolation  of  the  linear  relationship  found  for  the  lower  voltages. 

The  Tektronix  high  voltage  probe  with  oscilloscope  was  connected  to  one  of  the 
ceramic  insulator  conductors  in  order  to  measure  several  pulsed  voltage  characteristics.  For 
pulsed  voltages  in  the  range  of  12  to  30  kV,  the  voltage  decay  time  (one  measure  of  the  pulse 
width)  was  observed  to  range  from  1.385  to  1.600  ps.  A  typical  voltage  waveform  is  shown 
in  Figure  7,  and  the  plot  illustrates  the  typical  "ringing"  oscillatory  pattern  of  the  voltage 
decay.  The  resonance  period  for  these  oscillations  was  found  to  be  100  ns  and  this  was 
independent  of  the  applied  voltage.  These  oscillations  are  typically  observed  in  the  literature 
(Mizuno  et  al..  1993;  Masuda  and  Nakao,  1990;  Creyghton,  1994). 

Pulsed  corona  onset  was  found  to  occur  between  an  input  voltage  of  13  and  14  kV. 
The  corresponding  peak  voltages  seen  on  the  oscillograms  were  1.04  and  23.4  kV, 
respectively.  Below  onset,  the  output  voltage  reading  was  very  low  and  had  erratic  and  rapid 
variations.  Once  onset  occurred,  the  output  voltage  rapidly  rose  to  23.4  kV.  Between  23.4 
and  41.0  kV  peak  pulsed  voltages  the  rise  time  and  decay  time  were  seen  to  vary  from  45  to 
60  ns  and  1.535  to  1.665  ps,  respectively.  It  was  found  that  there  was  no  significant 
difference  between  the  voltage  required  to  initiate  the  corona  and  that  required  to  maintain  it, 
i.e.,  there  was  no  hysteresis. 

Analysis  of  the  voltage  waveform  at  the  HV  line  side  of  the  capacitor  during  repetitive 
spark  gap  closures  was  performed  to  determine  the  voltage  loss.  It  was  observed  that  the 
frequency  of  the  capacitor  discharge  matched  the  frequency  of  the  rectified  wave  (60  Hz)  as 
expected,  and  the  discharge  was  very  fast,  1.6  ps,  compared  to  the  calculated  RC  capacitor 
charging  time  of  5  ms.  However,  there  was  a  loss  of  about  1.2  -  2.0  kV  once  the  capacitor 
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reached  its  maximum  charge.  This  loss  has  been  attributed  to  a  slow  leak  to  ground  in  the 
capacitors.  From  the  waveform,  it  was  observed  that  the  oscillations  present  in  the  output 
voltage  waveform  are  not  present  at  the  line  side  of  the  capacitor.  This  suggests  that  the 
source  of  the  oscillations  is  probably  due  to  inductance  in  the  circuit. 

Load  resistors  (R2  in  Figure  2)  are  present  in  the  circuit  since  the  impedance  of  the 
wire-cylinder  separation  gap  is  too  high  to  sustain  a  pulsed  discharge  in  the  gas  phase  reactor. 
Generally,  the  value  of  this  resistance  is  100  -  200  ohms.  The  number  of  load  resistors  (8 
ohms  each)  were  changed  to  analyze  the  effect  of  the  total  resistance  on  the  output  voltage 
waveform.  The  decay  time  was  strongly  affected  by  the  load  resistance.  All  subsequent  NO 
removal  experiments  were  performed  with  approximately  150  ohms  tail  resistance  in  place. 

3.2.  Pulse  Waveform  Characteristics 

Figures  7,  8,  and  9  show  typical  pulsed  voltage,  current,  and  power  characteristics 
measured  at  the  reactor  discharge  conductor.  Measurements  were  repeated  for  input  voltage 
dial  settings  of  30  kV  to  50  kV.  The  oscilloscope  was  set  to  average  the  waveform  over  100 
pulses.  The  amount  of  oscillation  in  the  voltage  waveform  was  observed  to  decrease  as  the 
voltage  was  increased.  The  current  pulses  are  generally  shorter  than  the  voltage  pulses,  and 
the  current  waveforms  also  show  a  reversal  to  negative  current  for  a  short  period  following  the 
peak  in  the  pulse.  Multiplying  a  corresponding  voltage  and  current  waveform  gives  a  power 
waveform;  the  power  waveforms  follow  a  trend  similar  to  that  of  the  current  pulses,  however, 
they  generally  had  a  double  peak. 

Table  I  summarizes  the  data  for  the  power  supply  characterization.  The  voltage  and 
current  were  measured  at  the  reactor  HV  electrode  using  the  probes  connected  to  the  digital 
storage  oscilloscope.  The  dial  setting  on  the  power  supply  was  varied  from  30  to  50  kV.  The 
peak  pulsed  voltages  are  also  shown  in  this  table  and  generally  increase  with  increasing  input 
power.  The  pulse  width  measured  at  1/2  maximum  for  the  voltage  pulse  varied  from  450  ns 
to  694  ns.  The  peak  current  was  also  seen  to  generally  increase  with  input  voltage,  although 
the  width  of  the  current  pulse  at  the  baseline  was  seen  to  be  invariant  with  the  applied  voltage. 

The  energy  per  pulse  was  calculated  by  integrating  the  voltage  times  the  current  over  a 
single  pulse.  The  power  per  pulse  varied  from  32.1  mJ/pulse  to  107  mJ/pulse,  for  input 
voltages  of  30  kV  and  50  kV,  respectively.  The  power  delivered  to  the  reactor  was 
determined  by  multiplying  the  energy  per  pulse  by  the  number  of  pulses  per  second  (60  Hz). 
The  power  was  found  to  vary  from  1.9  to  6.4  watts,  for  voltage  inputs  of  30  to  50  kV.  The 
energy  delivered  per  pulse  was  somewhat  higher  than  that  reported  by  Tas  (1995)  (6  -7 
mJ/pulse),  however,  the  total  power  delivered  was  slightly  higher  than  that  reported  by  Tas 
(1995)  (0.05  to  1.5  watts).  The  energy  per  pulse  was  similar  to  that  reported  by  Mizuno  et  al. 
(1993)  (20  to  56  mJ/pulse),  and  the  power  delivered  was  slightly  smaller  than  the  discharge 
power  reported  by  Mizuno  et  al.  (1993)  (5  to  14  watts).  Mizuno  et  al.  (1993)  measured  both 
the  input  power  and  the  discharge  power  and  they  found  that  the  discharge  power  was  linearly 
related  to  the  input  power;  the  discharge  power  was  20  to  30%  of  the  input  power.  Masuda 
and  Nakao  (1986)  reported  a  power  output  in  the  range  of  0.5  to  6.0  watts. 
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Tas  (1995)  varied  the  frequency  of  his  power  generator  (5  -  240  Hz)  and  held  the 
voltage  at  a  constant  level  (28  kV;  16  kV/cm).  Mizuno  et  al.  (1993)  fixed  the  frequency  at 
250  Hz  and  varied  the  voltage  from  6.5  to  14  kV  (8.5  -  19  kV/cm).  Masuda  and  Nakao  (1986) 
varied  their  corona  input  voltage  from  0  to  60  kV  (0-12  kV/cm).  In  the  present  study,  the 
frequency  was  held  constant  at  60  Hz  and  the  voltage  was  varied  from  30  -  50.  Figure  10 
shows  that  the  power  output  of  the  FSU  pulsed  power  supply  was  approximately  linear  with 
the  input  voltage. 


Table  1 .  Power  Supply  Characteristics 
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It  is  important  to  note  that  the  operating  procedures  for  the  experiments  reported 
here  required  cleaning  the  stainless  steel  spark  gap  electrodes  after  each  run  in  order  to 
ensure  reproducibility  of  the  data.  A  light  sanding  with  fine  sandpaper  and  wiping  with 
ethanol  (see  Section  2.5)  was  performed  after  each  run  because  the  current  and  voltage 
characteristics  can  be  affected  by  buildup  of  deposits  on  the  electrodes.  Future  work  will 
consider  running  the  spark  gap  in  an  inert  or  nitrogen  gas  in  order  to  minimize  this  buildup 
on  the  electrodes  and  to  reduce  the  production  of  ozone  by  the  spark  gap.  It  can  be  noted 
that  Mizuno  et  al.  (1993)  used  tungsten-coated  electrodes. 

3.3.  Pulsed  Corona  Reactor  Experiments  with  Dry  Air  and  NO 

Figure  1 1  shows  the  steady-state  NO  concentration  vs.  applied  voltage  for  the  three 
different  residence  times  used  in  this  study  (Run  #'s  3,  4,  13,  14,  22,  23,  10,  11,  20,  21,  16,  17, 
5,  7,  26,  27,  18,  19).  The  gas  composition  was  dry  air,  150  ppm  NO,  and  no  H20  or  ethylene. 
It  is  clear  from  this  figure  that  at  30  kV  the  removal  of  NO  is  very  low  at  the  short  residence 
times.  There  is  also  a  sharp  transition  between  30  kV  and  40  kV  at  these  short  residence  times 
and  this  behavior  is  consistent  with  that  observed  by  Mizuno  et  al.  (1993)  and  Masuda  and 
Nakao  (1990).  At  a  residence  time  of  0.6  minutes,  the  maximum  NO  removal  (@  50  kV)  was 
75%. 
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Table  2,  Summary  of  Experimental  Results 


Table  2.  Summary  of  Experimental  Results  (Continued) 


Table  2.  Summary  of  Experimental  Results  (Concluded) 
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Figure  12  shows  the  NO  concentration  vs.  residence  time  in  the  reactor  for  various 
voltage  dial  settings  (Run  #'s  3,  4,  13,  14,  22,  23,  10,  11,  20,  21,  16,  17,  5,  7,  26,  27,  18,  19). 
The  gas  composition  was  the  same  as  described  above.  It  is  again  clear  that  at  the  low 
voltage,  i.e.,  30  kV,  the  removal  of  NO  at  short  residence  times  is  very  small  as  mentioned 
above.  At  the  higher  voltages  the  removal  is  much  greater.  At  the  long  residence  time  of  12 
min.  virtually  100%  of  the  NO  was  removed  for  the  two  higher  voltages  of  43.9  and  48.8  kV. 
It  is  important  to  note  that  at  this  long  residence  time  the  N02  also  decreased.  As  the  flow  rate 
is  increased  and  the  residence  time  is  decreased,  a  corresponding  drop  in  the  amount  of  NO 
removed  is  observed.  It  is  also  important  to  note  that  the  efficiency  of  NO  removal  is 
improved  at  the  higher  flow  rate  as  shown  in  Table  2.  Here  the  amount  of  energy  to  remove 
one  mole  of  NO  is  approximately  100  eV  for  the  higher  voltages  (42  to  60  kV)  at  10,000 
seem.  This  number  can  be  compared  to  that  obtained  by  Tas  (1995).  For  dry  nitrogen  with 
6%  oxygen  and  1200  ppm  NO,  Tas  obtained  an  energy  of  approximately  850  eV/mole  of  NO 
removed.  He  did  not  report  data  at  the  higher  oxygen  content  used  in  the  present  study, 
however  he  did  observe  an  increase  in  energy  per  NO  as  oxygen  content  increased  from  0  to 
6%. 


Figure  12  also  shows  an  approximately  linear  decrease  in  NO  concentration  with 
residence  time  for  the  lower  voltage  condition  (30  kV,  Run  #'s  5,  7,  26,  27,  18,  19).  This  is 
consistent  with  an  approximate  relationship  derived  by  Tas  (1995).  At  higher  power  input, 
with  low  NO  input,  Tas  (1995)  found  a  first  order  decay  in  NO  concentration,  and  this  trend  is 
also  seen  approximately  in  the  data  of  Figure  12  for  the  higher  voltages  (40  and  50  kV). 

Analysis  of  the  effluent  from  the  reactor  using  GC/MS  indicated  that  about  15  ppm 
N20  was  produced  in  all  the  cases  with  NO  feed  in  dry  air.  This  approximate  amount  of  N20 
appears  to  be  formed  in  the  system  regardless  of  the  presence  of  other  species,  i.e.,  water  or 
ethylene.  The  independence  of  N20production  on  feed  composition  is  consistent  with 
observations  reported  by  Mizuno  et  al.  (I  993). 

3.4  Pulsed  Corona  Reactor  Experiments  with  Dry  Air,  NO,  and  Ethylene 

Figure  13  shows  the  NO  concentration  at  the  reactor  outlet  as  a  function  of  the  applied 
voltage  for  the  case  with  500  ppm  ethylene  and  for  the  case  without  any  ethylene  (Run  #'s  18, 
19,  28,  29,  13,  14,  16,  17,  33,  34,  36,  37,  49,  50,  51,  52,  53,  54,  55,  56,  57,  58,  59,  60).  It  can 
be  seen  that  at  lower  voltages  (30  kV)  the  difference  between  the  two  cases  is  small.  The 
largest  difference  is  seen  at  40  kV  where  the  presence  of  ethylene  causes  the  NO  removal  to 
decrease  from  approximately  50%  to  approximately  80%.  At  the  highest  voltage,  50  kV,  the 
presence  of  ethylene  causes  the  NO  removal  to  go  from  about  80%  to  about  90%. 

Energy  consumption  calculations  in  the  present  report  use  the  relationship  (Tas,  1995) 

1  eV  /  molecule  NO  ~  96.4  kJ  /  mole  NO  removed.  (37) 

The  energy  consumption  per  NO  removed  at  35,  40,  45,  and  50  kV  was  approximately  50 
eV/NO,  and  at  30  kV  it  was  approximately  250  eV/NO.  Mizuno  et  al.  (1993)  found  the 
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energy  consumption  to  vary  from  80  to  151  eV/NO  removed  when  ethylene  was  added  to  the 
reactor  for  a  dry  gas  with  10%  oxygen  and  200  to  800  ppm  NO.  The  lower  value  of  80 
EVNO  compares  well  to  the  value  of  50  eV/NO  obtained  in  the  present  study. 

Analysis  by  GC/MS  of  the  effluent  from  the  reactor  with  ethylene  in  the  gas  feed 
indicates  that  at  45  kV,  23%  of  the  ethylene  has  reacted  to  products  and  that  at  50  kV,  26%  of 
the  ethylene  has  reacted.  This  observation  is  based  on  a  decrease  in  the  ethylene  peak  area 
measured  by  the  GC/MS.  Since  the  feed  gas  contains  500  ppm  ethylene  this  indicates  that 
about  115  to  130  ppm  ethylene  is  decomposed  by  the  pulsed  corona.  These  measurements 
were  made  for  the  reactor  run  with  a  0.6  min  residence  time  and  with  170  to  200  ppm  NO 
(Run  #s  58,  59,  60).  The  corresponding  amounts  of  NO  removed  from  the  reactor  under  these 
conditions  were  140  and  170  ppm.  Therefore,  a  slightly  larger  amount  of  NO  was  removed 
than  would  be  expected  if  the  reaction  had  a  1:1  stoichiometry.  No  other  reaction  products 
were  detected  with  the  GC/MS,  although  additional  experimentation  with  other  sampling 
techniques  and  analysis  may  be  necessary  to  validate  this  observation. 

3.5.  Pulsed  Corona  Reactor  Experiments  with  Humid  Air,  NO,  and  Ethylene 

Table  3  compares  the  NO  removal  for  the  four  cases  of  NO  only,  NO  with  ethylene, 
NO  with  water,  and  NO  with  water  and  ethylene.  All  runs  were  made  at  10,000  seem  and  50 
kV  voltage  dial  setting.  The  concentration  of  NO  in  the  feed  was  in  the  range  of  140  to  180 
ppm.  The  gas  was  at  room  temperature  (approximately  25°  C),  and  the  water  content  was  no 
higher  than  1  -  2%  by  weight  (saturation  at  room  temperature,  or  3  mole  percent).  Without 
additives  the  NO  removal  was  about  75%  (average  of  Run  #'s  16,  17,  39,  40).  The  effect  of 
500  ppm  ethylene  was  more  dramatic  than  the  effect  of  water,  i.e.,  93%  removal  (average  of 
Run  #s  36,  37,  49,  50,  60)  vs.  81%  removal  (average  of  Run  #'s  47,  48).  In  addition,  the  case 
with  water  and  ethylene  showed  only  a  relatively  small  percentage  improvement  over  the  case 
with  ethylene  only,  i.e.,  98%  removal  (average  of  Run  #'s  41,  42,  43,  45)  vs.  93%  removal. 

Table  3  also  shows  the  energy  costs  for  removing  one  molecule  of  NO  for  each  of  the 
four  conditions  reported.  The  NO  only  case  required  76  eV/NO  to  remove  73%  of  the  input 
NO.  In  dry  nitrogen  containing  10%  oxygen  the  data  of  Mizuno  et  al.  (1993)  indicate  290  to 
370  eV/NO,  and  in  dry  nitrogen  containing  0  to  6%  oxygen  the  data  of  Tas  (1995)  indicate 
210  to  850  eV/NO.  With  the  addition  of  ethylene  the  present  study  suggests  56  eV/NO  while 
the  data  of  Mizuno  et  al.  (1993)  for  10%  oxygen  in  nitrogen  gives  80  to  150  eV/NO.  For  the 
case  of  10%  oxygen  in  nitrogen  with  either  4.5  or  9.9%  water,  Mizuno's  data  gives  46  to  74 
eV/NO,  independent  of  water  content.  These  latter  values  are  in  the  general  range  of  those 
given  in  Table  3. 
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Table  3.  Effect  of  Water  on.  NO  Removal  (50  kV,  10,000  scan  air) 


NO  only 

NO  +  ethylene 

NO  +.  water 

NO  +  ethylene 
+  water 

75% 

93% 

81% 

98% 

76  eV/NO 

56  eV/NO 

66  eV/NO 

58  eV/NO 

For  all  experiments  run  with  the  feed  air  bubbled  through  a  water-filled  chamber,  it 
was  observed  that  a  slight  film  was  formed  on  the  inside  of  the  glass  sampling  tubes  when  gas 
samples  were  taken.  This  film  was  dissolved  in  water  and  the  pH  of  the  resulting  solution  was 
measured.  No  decrease  in  pH  was  found,  although  again  further  validation  of  this  finding 
needs  to  be  made.  Further  work  needs  to  be  performed  to  identify  the  composition  of  this 
material  and  to  detect  other  reaction  products  from  the  corona  reactor. 

In  order  to  compare  data  obtained  under  a  variety  of  conditions  in  our  experiments  to 
data  from  several  literature  sources.  Figures  14  and  15  show  "percentage  NO  removal"  for 
low  initial  concentration  cases  and  "ppm  NO  removed"  for  high  concentration  cases, 
respectively,  as  functions  of  a  factor  defined  as  "the  power  per  volume  (energy  input  divided 
by  flow  rate)  divided  by  the  square  root  of  the  residence  time".  This  type  of  plot  was 
suggested  by  Masuda  and  Nakao  (1986),  where  they  were  able  to  compress  data  from  a  range 
of  power  and  residence  time  experiments  onto  a  single  plot.  No  theoretical  justification  for 
this  plot  is  available  at  present,  however  it  serves  as  a  convenient  means  to  compare  data 
obtained  under  a  variety  of  conditions. 

Figure  14  shows  that  the  data  from  the  present  study  for  dry  air  falls  along  a  single 
smooth  curve.  There  is  one  outlier  point  not  shown  in  the  figure  that  arises  at  low  voltage  (30 
kV)  and  short  residence  time;  this  may  be  due  to  the  low  voltage  being  insufficient  to  initiate 
the  expected  chemical  reactions.  Data  from  Masuda  and  Nakao  (1986)  and  Ohkubo  et  al. 
(1994),  as  well  as  data  from  the  present  study  for  the  cases  with  ethylene  present  and  water 
present  are  shown  for  comparison.  Masuda  and  Nakao  (1986)  considered  the  case  of  dry  air 
with  180  ppm  NO  and  180  ppm  NH3.  The  percentage  NO  removal  for  this  case  with 
ammonia  is  higher  than  that  of  our  data  for  dry  air.  Our  data  for  the  cases  of  ethylene  and 
water  compare  well  to  the  data  of  Masuda  and  Nakao  (1986)  at  1.8  Wh/m3-s1/2.  The  data  of 
Ohkubo  et  al.  (1994)  for  cases  of  20%  02  in  nitrogen  fall  somewhat  below  our  data  and  their 
data  for  the  case  with  ammonia  does  not  show  much  of  an  effect  of  ammonia  on  NO  removal. 
It  can  be  concluded  from  this  figure  that  the  general  range  of  performance  of  the  reactor 
developed  in  this  project  is  similar  to  other  reactors  reported  in  the  literature.  Because  of  the 
wide  range  of  gas  feed  conditions  and  possibly  variation  in  methods  of  determining  power 
output,  precise  comparison  between  different  studies  is  difficult. 

Figure  15  shows  data  from  Tas  (1995)  and  Mizuno  et  al.  (1993)  for  much  higher  levels 
of  initial  NO  concentration.  The  data  obtained  by  Mizuno  et  al.  (1993)  for  10%  O2  with  4.5% 
or  9.9%  H20  in  nitrogen  is  seen  in  this  figure  to  be  very  close  to  that  obtained  by  Tas  (1995) 
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for  5%  02  with  11%  H20  in  nitrogen.  It  is  interesting  to  note  that  the  data  of  Mizuno  et  al. 
(1993)  shows  no  effect  on  reactor  performance  of  water  content  in  the  range  of  4.5  %  to  9.9%. 
The  data  of  Mizuno  et  al.  (1993)  also  shows  a  large  effect  of  ethylene  on  NO  removal  from 
nitrogen  gas  with  10%  02.  Ethylene  removal  data  (Mizuno  et  al.,  1993)  shows  a  significantly 
larger  energy  consumption  than  the  work  reported  here.  It  is  not  clear  why  that  shady  differs 
from  the  present  work,  although  there  are  differences  in  the  gas  feed  composition  (oxygen 
content)  that  warrant  further  investigation.  The  "no  ethylene"  data  of  Mizuno  et  al.  (1993) 
also  show  significantly  higher  energy  costs  than  that  of  the  present  study. 
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4.  MODELING  RESULTS 


Although  prediction  of  the  performance  of  pulsed  streamer  corona  reactors  through 
mathematical  modeling  has  been  pursued  by  a  number  of  investigators,  no  theoretical  model 
has  been  tested  and  quantitatively  compared  with  actual  experimentally  measured  output  from 
a  pulsed  corona  reactor.  Empirical  models  developed  by  Masuda  and  Nakao  (1990)  describe 
the  removal  efficiency  as  a  function  of  the  input  power,  but  do  not  explicitly  account  for  the 
complex  chemistry  occurring  in  the  system.  Many  other  theoretical  models  have  been  applied 
to  corona  systems,  however,  only  qualitative  trends  have  been  predicted  due  to  the  nature  of 
the  assumptions  made  concerning  the  details  of  the  streamer  structure.  A  range  of  approaches 
have  been  made  to  modeling  N0X/S02  removal  (Wren,  1989;  Alekseev,  1993;  Civitano  et  al., 
1988;  Creyghton,  1994;  Dinelli  et  al.,  1990;  McFarlane  and  Wren,  1991;  Mizuno  et  al.,  1993; 
Tokunaga  and  Suzuki,  1990;  Gallimberti,  1988;  Masuda  and  Wu,  1987;  Masuda  et  al.,  1987; 
Masuda  et  al.,  1991)  and  ozone  synthesis  (Peyrous  et  al.,  1989). 

Models  of  pulsed  streamer  corona  can  be  classified  into  those  that  consider  a)  streamer 
propagation  through  species  continuity  equations  for  electrons  and  ions,  and  b)  batch  chemical 
reaction  kinetic  models.  The  analysis  of  streamer  propagation  is  generally  performed  with 
one  dimensional  convective/diffusion/reaction  models  that  consider  electron  propagation  in  a 
non-uniform  electrical  field  (see  Creyghton,  1994,  for  good  introductory  review).  One  of  the 
objectives  of  these  models  is  to  determine  the  minimum  voltage  required  for  streamer 
inception.  Batch  chemical  kinetics  models  usually  consider  a  large  number  of  reacting 
species,  including  electrons,  ions,  radicals,  and  molecular  species;  however,  the  electron 
concentration  is  input  either  as  a)  an  empirically  determined  function  (as  in  measured  current 
pulses),  or  b)  through  solution  of  the  Boltzmann  equation  for  electron  distribution  as  a 
function  of  the  applied  electrical  field.  Coupling  the  streamer  models  to  the  batch  kinetic 
models  has  been  attempted,  however,  the  approaches  that  have  been  used  require  severe 
limiting  assumptions  on  the  nature  of  the  streamer  structure  (i.e.,  ideal  uniform  disk 
structures).  It  is  therefore  clear  that  it  is  necessary  to  further  develop  the  stream er/chemistry 
models  in  order  to  have  solid  quantitative  comparison  of  experimentally  measured  reactor 
output  with  predictive  models. 

The  modeling  results  in  the  present  study  have  focused  on  the  use  of  a  basic  set  of 
chemical  reactions  shown  in  Table  4.  The  reaction  rate  constants  for  all  of  these  reactions, 
except  the  first  two  electron  reactions  with  N2  and  02,  were  taken  from  the  literature 
(McFarlane  and  Wren,  1991).  The  reaction  rate  constants  for  these  two  reactions  were  varied 
to  compare  with  the  experimental  data.  These  reactions  are  expected  to  be  highly  dependent 
upon  the  electric  field  strength.  Figures  16  and  17  show  the  NO  and  N02  as  a  function  of 
residence  time  for  various  values  of  this  initiation  reaction  constant.  It  can  be  seen  in  Figure 
16  that  at  the  lower  values  of  the  reaction  rate  constant  the  decrease  in  NO  with  residence  time 
is  much  slower  than  that  at  the  higher  values.  The  data  from  the  experiments  reported  in  the 
present  work  showed  a  NO  removal  at  36  seconds  of  about  74%.  This  would  fall  between  the 
profiles  computed  with  k,  =  2  x  1012  cm3/mole-s  and  kj  =  5  x  1012  cm3  /mole-s  as  shown  in  this 
figure.  Figure  17  shows  the  corresponding  N02  profiles.  It  is  very  interesting  to  note  that  at 
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Table  4.  Reactions  Used  to  Compute  the  Different  Concentration  Profiles 

ELEMENTS  H  E  O  N  C  END 

SPECIES  K2  H  02  O  0*  OH  H02  H202  H20  N  N2  NO  03  HN02  HN03  E  N02  NC>3  N20 
N204  C2H4  HCHO  HCOOH  CO  C02  02H3O  HOC2H4  B0C2H402  H0C2H40  CH20H  HOCH2CHO 
02CH2OH  0CH20H  HCO  END 


REACTION 

REACTION  RATE  CONSTANT 

N2  +  E  =>  2N  +  E 

0./SB+13 

02  +  E «>  0  +  0*  +  £ 

0.5E+13 

mo  +  E  =>  H  +  OH  +  E 

0,5 E+13 

N  +  N  N2 

6 . 384E+10 

N  +  0  =>  NO 

1  +  445B+11 

N  +  02  =>  KO  +  O 

53604700 

Ml — — 

2..047E+13 

N  +  NG2  2N0 

1 . 8069E+12 

N  +  N02  ->  N2  +  02 

1 . 0841E‘Tl3 

N  +  NC2  N2  +  20 

1.385E+12 

N  +  OH  =>  NO  +  H 

- -1.4458+13 

N  +  N02  »>  N20  +  6 

1 . S069E+12 

0*  *>  0 

6.95E+8 

0*  +03  =>  202 

7.2275E+13 

0*  +  03  20  +  02 

7.2276E+13 

O*  +  N20  N2  +  02 

2 . 9513E+13 

0*  +  N20  2N0 

4 . 03 54 S+ 13 

0*  +  H2  H  +  OH 

6 . 023E+13 

0*  +  H20  2 OH 

1 . 3251E+14 

0  +  6  =>  02 

1.0239E+11 

O  +  H  ->  OH 

1 . 0B4I4E+12 

b  4-  02  =e>  03 

7 » 829913+9 

0  +  03  202 

:  4 . S184E+9 

0  +  H20  ->  OH  +  02 

3.433E+13 

0  +  OH  =>02  +  H 

1,  987E+13 

OH  +  H2  «  H20  +  H 

4.035E+3 

OH  +  OH  =>  H20  +  0 

1 . 0841E+12 

OH  +  OH  =>  H202 

9 937E+12 

OH  +  H02  »>  H20  +  02 

6.62S3E+13 

OH  +  H202  =>  H20  +  H02 

X .  0239E+12 

OH  +  H  H20 

1.4452E+X3 

H  +  02  =>  H02 

8. 4928+11 

H  +  03  ->  OH  +  02 

1 . 6S64E+13 

K  4  HD2  *>  H2  +  02 

3 . 372E+I2 

H  +  H02  «>  2 OH 

4.3365E+13 

H  +  H02  =>  H20  +  6 

1.44552E+12 

H  +  H202  «>  H20  +  OH 

6 . 023E+10 

H02  +  H02  ->  H202  +  02 

1 , 0239E+12 
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NO  +  OH  «>  HN02 


NO  +  03  «>  N02  +  02 


NO  +  O  >  N  +  02 


NO  +  O  N02 


NO  +  NO 3  2N02 


N02  +  OH  =>  HN03 
NO 2  +  03  »>“N03  +  02 


N02  +  O  NO  +  02 


NO 2  +  OH  =>  HO 2  +  NO 
HN02  +  OH  N02  +  H20 


HN03  +  OH  .>  H20  +  N03 


9.0345E+12 
"s .4322E+9 


2 . 77058E+9 


1 ,  3973E+12 


■1-7466E+13 


.1 -4455E+X3 
*  192736  00 


5-;.8423E+l2 


1.38529B+13 
2 . 951E+12 


9  *  034E+10 


H204  4-  M  ->  2N02  +  M 


2N03  =>  2N02  +  02 


N03  +  O  »>  N02  +  02 


2N02  +  M  ->  N204  +  M 


N02  4*  N03  «>  N02  +  NO  +  02 


2.0639E+9. 
1  9177E+8 


X-  0239E+13 


NO  4-  N03  =>  2ND  +  02 

1 , 8352E+11 

03  4  N  =>  NO  +  02 

6 , 023E+7 

3 . 0553E+14 
"£  +  4529B+8 


OH  +  H  ->  K2  +  O 
H20  4-  O” =>20H 


1.964E+7 


I , 1249E+12 


H202  +  H  H02  +  H2 

3 . 121.9S+9 

H202  +  0  =>  OH  +  H02 

_  : 

1 . 0262E+9 

~>  OH  +  202 


03  +  OH  ->  H02  +  02 


HN02  +  0  =>  NO 2  +  OH 


HN03  +  HN02  =>  2N02  4-  H20 


HN03  +  O  *>  N03  +  OK 


HN03  +  H  »>  HN02  +  OH 


N03  +  H  =>  N02  +  OH 


N03  *  OH  =>  N02  +  HQ 2 


N03  +  H02  «>  HN03  +  02 


NOB  +  H02  =>  N02  +  OH  +  02 


» >  NO  +  OH 


O  +  C2H4  «>  C2H30  +  H 


«>  HCHO  +  CO  4*  OH 


C2K4  4-  OH  **>  BOC2H4 


HOC2H4  +  02  =  >  H0C2H402 


HOC2H402  4*  NO  «>  H0G2H40  +  N02 


H0C2H40  HCHO  +  CH20H 


H0C2H4O  +  02  «>  H0CH2CH0  +  H02 


CH20H  +  02  =>  HCHO  +  H02 
H02  4-  HCHO  02CH20H 


02CH20H  +  NO  »>  0CH20H  4-  N02 


0CH20H  +  02  =>  HCOOH  +  H02 


1 ♦ 1068E+9 
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the  highest  k,  there  is  a  maximum  in  N02  concentration,  and  that  at  longer  residence  times  the 
N02  decreases.  This  maximum  in  N02  is  observed  in  the  data  of  Masuda  and  Nakao  (1990) 
with  600  ppm  of  ammonia  present. 

For  the  highest  value  of  k,  used  in  Figures  16  and  17  (k,  =  5  x  1012  cm3/mole-s). 
Figure  18  shows  the  NO,  N02,  and  N20  profiles.  This  figure  indicates  that  the  increase  in 
N20  is  very  small  until  the  residence  time  is  long,  and  that  all  of  the  removal  of  NO  cannot  be 
accounted  for  in  the  formation  of  N02  and  N20;  i.e.,  the  sum  of  N02,  N20  and  all  nitrogen 
oxide  species  does  not  equal  the  initial  NO  gas  feed  to  the  reactor.  This  indicates  that  the 
conversion  of  some  N02  to  N2ando2may  account  for  the  nitrogen  balance. 

Figures  19,  20,  and  21  show  the  corresponding  profiles  (NO,  N02,  and  HN03)  for  the 
cases  with  1%  water  present  and  the  same  variation  in  rate  constants  as  used  in  the  previous 
"no  water"  results.  The  trends  are  similar  to  those  shown  in  the  previous  figures,  however,  the 
formation  of  HN03  can  now  be  observed.  Comparison  of  Figure  16  with  Figure  19  shows 
that  the  presence  of  1%  water  vapor  significantly  enhances  the  NO  removal.  For  example,  the 
case  of  k,  =  5  x  10n  shows  a  22%  removal  for  the  dry  case  and  a  36%  removal  for  the  wet 
case.  Comparison  of  Figure  17  with  Figure  20  shows  a  shift  in  the  time  where  the  NO 
maximum  occurs.  For  example,  in  the  case  of  k3  =  5  x  1012,  the  NO  maximum  occurs  at  about 
30  seconds  in  the  dry  system  and  the  maximum  occurs  at  about  18  seconds,  Figure  20,  in  the 
wet  system.  From  comparison  of  Figure  17  and  Figure  22,  it  is  also  obvious  that  the 
concentration  of  N02  is  much  lower  in  the  wet  system.  Figure  23  shows  the  increase  in  HN03 
and  again,  the  quantitative  amount  of  products  formed  as  the  NO  decomposes  also  indicates 
that  some  nitrogen  from  NO  may  be  converted  to  N2. 

The  presence  of  500  ppm  ethylene  in  the  dry  system  (Figure  24)  gives  about  99% 
removal  of  NO  at  about  50  seconds  residence  time.  This  can  be  compared  to  the  greater  than 
60  second  residence  time  required  in  the  no  additive  case  shown  in  Figure  18.  The  shape  of 
the  N20  and  03  curves  in  Figure  25  are  also  slightly  different  than  those  for  the  no  additive 
case.  Figure  26  shows  that  500  ppm  ethylene  strongly  enhances  the  production  of  ozone. 

Figure  27  summarizes  the  four  cases  of  dry  air,  air  with  ethylene,  air  with  water,  and 
air  with  ethylene  and  water.  With  the  rate  constants  used  in  the  simulation,  all  four  cases 
appear  to  reduce  NO  from  initial  concentrations  by  over  99%  in  60  seconds.  At  36  seconds 
the  dry  air  case  showed  a  90%  reduction  in  NO,  the  air  with  ethylene  case  showed  a  95% 
reduction  in  NO,  the  air  with  water  case  showed  a  98%  reduction  in  NO,  and  the  air  with 
water  and  ethylene  case  showed  a  100%  reduction  in  NO.  The  slowest  rate  of  reduction  is 
that  for  the  dry  air  alone,  as  observed  in  the  experimental  results,  although  the  experimental 
data  still  indicated  about  a  74%  drop  in  NO  with  this  residence  time.  The  next  slowest  case 
shown  in  the  model  is  for  the  case  of  air  with  ethylene.  This  does  not  correspond  to  the 
experimental  results  where  the  air  with  water  case  showed  a  slower  rate  of  NO  removal  than 
the  air  with  ethylene  case  —  81%  verses  93%,  respectively,  at  a  36  second  residence  time. 
This  discrepancy  may  be  due  to  differences  in  the  rate  constants  used  in  the  model  vs.  those  in 
the  experimental  system,  or  to  possible  differences  in  the  amount  of  water  used  in  the 
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experimental  study.  A  further  problem  with  the  model  is  that  the  reaction  of  ethylene  with 
ozone  (according  to  the  literature)  is  five  orders  of  magnitude  slower  than  that  of  ethylene 
with  hydroxyl  radicals  or  singlet  oxygen.  In  addition,  the  exact  rate  constants  for  all  of  the 
individual  steps  in  the  ozone  reaction  mechanism  are  not  known.  The  computer 
simulation  therefore  does  not  include  the  ozone  reaction  with  ethylene.  The  simulation 
shows  that  the  case  of  air  with  water  and  ethylene  has  an  initially  slower  slope  than  the 
case  with  water  alone,  however,  these  curves  cross  at  about  30  seconds,  and  the  case  of  air 
with  water  and  ethylene  drops  to  zero  at  a  shorter  residence  time  than  the  case  of  air  with 
water. 


Figure  28  shows  that  the  addition  of  water  dramatically  affects  the  ethylene 
concentration,  and  that  without  water  very  little  ethylene  is  removed.  The  experimental 
results  show  that  ethylene  is  removed  significantly  (approximately  25%),  even  without 
water  present  in  the  system.  This  again  may  be  due  to  our  neglecting  the  ozone-ethylene 
reactions  as  discussed  above.  Figures  29  and  30  show  two  of  the  hydrocarbon  byproducts 
produced  from  ethylene  decomposition  in  the  pulsed  corona  discharge  for  the  cases  of  air 
with  and  without  water.  The  production  of  formaldehyde  and  glycolaldehyde  are  both 
strongly  enhanced  by  the  addition  of  water. 
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5.  CONCLUSIONS 


The  major  conclusions  from  this  study  include: 

1.  The  addition  of  ethylene  to  dry  air  in  the  pulsed  corona  reactor  results  in 
significant  enhancement  of  NO  removal  With  a  residence  time  of  0.6  min  and  a  power 
input  of  6.4  watts,  74%  of  the  NO  was  removed  from  dry  air  at  a  cost  of  76  eV/NO.  Upon  the 
addition  of  500  ppm  ethylene  under  the  same  reactor  conditions,  93%  of  the  NO  was  removed 
at  a  cost  of  56  eV/NO.  The  qualitative  effect  of  ethylene  confirms  data  reported  in  the 
literature  (Mizuno  et  al.,  1993).  The  data  of  Mizuno  et  al.  (1993),  however,  indicate  an  energy 
cost  of  between  80  and  150  eV/NO  with  the  presence  of  500  ppm  ethylene.  The  upper 
number  corresponds  to  200  ppm  NO  removed  and  the  lower  number  corresponds  to  800  ppm 
ethylene  removed.  Quantitative  analysis  by  GC/MS  indicates  that  roughly  equivalent 
amounts  of  ethylene  are  decomposed  for  NO  removal.  This  result  has  not  been  previously 
reported  since  the,  only  reported  work  on  the  effects  of  ethylene  (that  of  Mizuno  et  al.  1993) 
did  not  quantify  the  amount  of  ethylene  decomposed.  It  may  however,  be  inferred  from  the 
data  of  Mizuno  et  al.  (1993)  that  a  1:1  correspondence  is  not  to  be  expected  since  they  were 
able  to  remove  about  800  ppm  NO  with  an  addition  of  500  ppm  ethylene. 

2.  The  presence  of  water  in  air  containing  ethylene  results  in  a  small  increase  in  NO 
conversion  efficiency.  About  93%  of  the  NO  was  removed  at  a  cost  of  56  eV/NO  for  the 
ethylene  only  case,  as  mentioned  above.  Upon  the  combination  of  water  and  ethylene,  98%  of 
the  NO  was  removed  at  a  cost  of  58  eV/NO.  Water  presence  in  air  containing  NO  only  (no 
ethylene)  gives  an  NO  removal  of  81%  at  a  cost  of  66  eV/NO.  The  combined  effects  of 
ethylene  and  water  vapor  have  not  been  previously  reported  in  the  literature.  Analysis  by 
GC/MS  was  not  able  to  detect  other  ethylene  reaction  products  that  would  be  expected  based 
upon  known  chemistry.  Mizuno  et  al.  (1993)  detected  acetic  acid  reaction  products  using 
FTIR  spectroscopy.  The  reaction  pathways  of  known  atmospheric  chemistry  do  not  appear  to 
directly  favor  the  production  of  acetic  acid  as  shown  in  the  model  simulation  results  and  the 
basic  chemistry  obtained  from  atmospheric  chemistry  literature;  this  work  tends  to  show  more 
formic  acid  and  formaldehyde  produced.  Further  work  is  recommended  to  search  for 
formaldehyde  (GC/MS)  and  acetic  and  formic  acids  (HPLC).  The  model  results  indicate  that 
water  strongly  enhances  the  formation  of  acetic  acid  and  glycolaldehyde. 

3.  The  removal  of  NO  from  dry  air  shows  a  sharp  transition  above  30  kV  for  the  shorter 
residence  times  reported  in  this  work.  This  sharp  transition  in  removal  efficiency  has  been 
reported  by  Mizuno  et  al.  (1994),  however  the  quantitative  value  of  the  power  input  in  the 
present  system  differs  from  that  reported  by  Mizuno  et  al.  (1995). 

4.  Results  of  the  present  study  show  good  comparison  with  those  obtained  by  several 
other  investigators.  For  low  NO  removal  (0  to  200  ppm),  the  data  in  the  present  study  show  a 
similar  order  of  magnitude  in  energy  costs  to  that  of  Masuda  and  Nakao  (1986)  and  Ohkubo  et 
al.  (1994).  For  high  NO  removal  (200  to  1000  ppm),  the  data  of  Mizuno  et  al.  (1993)  show 
somewhat  higher  energy  costs  for  the  case  of  ethylene  present  in  comparison  to  our  data.  At 
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the  high  NO  removal  range,  the  data  of  Mizuno  et  al.  (1993)  compares  well  to  the  data  of  Tas 
(1995)  for  the  case  of  water  present. 

5.  Mathematical  simulation  of  the  chemical  reactions  occurring  for  NO  removal  indicates 

that  N02  initially  goes  up  at  the  shorter  residence  times,  reaches  a  maximum  at  an 
intermediate  residence  time,  and  thereafter  drops  at  the  long  residence  times.  The 
experimental  data  in  this  work  reflects  this  conclusion  because  at  the  long  residence  time  of 
12  minutes,  N02  was  found  to  decrease  to  very  low  levels  in  comparison  to  the  levels  found  at 
0.6  and  1.2  minutes.  The  model  also  indicates  that  ethylene  reduces  the  magnitude  of  the  N02 
peak,  however,  it  does  not  appear  to  shift  the  residence  time  where  the  peak  occurs.  The 
stoichiometry  of  the  reaction  products  from  the  model  simulation  indicates  that  a  small 
amount  of  NO  goes  to  produce  N20,  however,  at  long  residence  times  a  significant  quantity  of 
NO  may  react  to  form  N2.  GC/MS  data  confirms  the  presence  of  N20  in  the  trace  quantities 
(<  15  ppm)  predicted  by  the  model.  The  model  also  indicates  that  some  fraction  of  the  NO  is 
converted  to  N2  even  in  the  cases  with  water  and  ethylene  present. 

6'  The  effect  of  water  and  ethylene  on  NO  removal  in  the  model  simulations  show  similar 
trends  to  those  found  in  the  experimental  studies.  The  model  shouts  that  NO  is  least 
effectively  removed  in  air  alone,  and  that  NO  is  most  effectively  removed  in  air  with  1% 
water  and  500  ppm  ethylene.  The  model  also  shows,  however  a  more  pronounced  effect  of 
water  (at  1%)  than  of  ethylene  on  NO  removal,  while  the  experimental  study  shows  a  greater 
effect  on  NO  removal  from  the  presence  of  ethylene  than  of  water.  The  experiments  indicate 
that  air  with  water  showed  a  NO  removal  of  81%  and  that  air  with  ethylene  showed  an 
increased  NO  removal  of  93%.  The  model  simulations  indicate  that  1%  water  increases  the 
removal  to  98%  and  that  500  ppm  ethylene  increases  the  NO  removal  to  95%.  The  water 
content  of  the  experimental  study  is  not  precisely  known,  although  it  can  be  no  larger  than 
about  3%  mole  fraction  and  may  be  somewhat  less.  A  humidity  detector  that  has  been 
ordered  but  not  yet  received  can  be  used  to  measure  water  content  of  the  gas  more  accurately. 
The  magnitude  of  the  model  predictions  for  the  water  results  does  not  correspond  exactly  to 
the  experiments  due  to  uncertainties  in  the  electron  initiation  reaction  rates,  possible 
uncertainties  in  the  amount  of  water  added  during  the  experimental  studies,  and/or  possible 
differences  in  reaction  rate  constants  from  those  reported  in  the  atmospheric  chemistry 
literature. 

7.  GC/MS  analysis  shows  that  about  15  ppm  N20  is  produced  under  all  gas  feed 
conditions  of  this  study.  The  mathematical  simulation  shows  a  small  variation  in  N20 
production  (5-10  ppm)  for  the  cases  of  air  with  ethylene,  air  with  water,  or  air  with  ethylene 
and  water,  however,  this  variation  may  be  within  experimental  error  of  the  GC/MS  analysis. 
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6.  RECOMMENDATIONS  FOR  FUTURE  WORK 


The  major  recommendations  for  further  investigation  of  the  results  presented  in  this 
study  are: 

1 .  Perform  further  parametric  analysis  of  the  gas  feed  components  including  water,  CO, 
and  C02.  Due  to  time  and  equipment  delivery  constraints  the  present  study  was  unable  to 
perform  experiments  with  these  added  gases.  CO  is  expected  to  be  especially  important  in 
light  of  reactions  known  from  atmospheric  chemistry. 

2.  Perform  more  extensive  analysis  of  byproducts  formed  in  the  system.  GC/MS 
analysis  combined  with  HPLC  should  be  performed  to  search  for  quantitative  amounts  of  such 
reaction  products  as  acetic  acid,  formaldehyde,  and  formic  acid.  Sampling  methods  can  be 
improved  using  liquid  nitrogen  or  acid  bath  traps  to  collect  larger  quantities  of  sample. 
Efforts  should  be  made  to  try  to  determine  how  much,  if  any,  NO  is  converted  to  nitrogen 
through  reactions  in  the  model  simulation.  Quantification  of  HN03 production  should  also  be 
performed. 

3.  Comparison  of  the  model  predictions  with  experimental  results  for  direct 
measurements  of  N02,  03,  CO,  and  other  reaction  products.  The  present  study  did  not  provide 
quantitative  information  on  N02  due  to  problems  with  N02  calibration  gases.  Further  work 
with  N02  is  thus  essential.  The  production  of  03  was  not  studied  in  detail  due  to  problems 
with  the  ozone  monitor,  however,  further  work  to  test  the  model  predictions  should  be 
performed. 

4.  Measure  electron-N2  reaction  rate  constants  using  experiments  with  only  NO  in  N2. 
Compare  these  rate  constants  to  those  from  the  model  and  to  those  used  in  the  present  study. 
Measure  electron-02  reaction  rate  constants  for  ozone  production  in  experiments  containing 
only  02.  Compare  these  reaction  rate  constants  with  those  from  the  model. 

5.  Improve  reactor  power  calculations  by  determining  the  difference  between  delivered 
power  and  output  power,  i.e.,  determine  amount  of  power  lost  in  the  circuit  following  the 
approach  of  Mizuno  et  al.  (1993),  or  that  of  Tas  (1995). 

Additional  recommendations  for  consideration  include: 

1.  Install  temperature  control  apparatus  to  allow  reactor  operation  under  a  variety  of 
temperatures. 

2.  Develop  optimization  methods  to  find  optimal  operating  conditions  for  the  reactor  in 
terms  of  gas  feed  composition  and  power  variation. 

3.  Install  equipment  to  analyze  the  formation  of  aerosol  particles  in  the  outlet  of  the 
reactor. 
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Figure  2.  Circuit  Diagram  for  Pulsed  Power  Supply 
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Figure  4.  Gas-Phase  Pulsed  Corona  Reactor 
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Figure  8.  Current  Waveform  for  40  kV  Dial  at  Reactor  Outlet 
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Figure  9.  Integrated  Power  for  40  kV  Dial  at  Reactor  Outlet 
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Figure  10.  Reactor  Power  verses  Input  Voltage 
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input  Voltage 


Figure  1 1.  Effect  of  Voltage  on  NO  Removal  from  Dry  Air 
(Initial  Concentration  of  NO  150  ppm,  H,0  Absent,  Ethylene  Absent) 
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Figure  13.  Effect  of  Voltage  on  NO  Removal  With  and  Without  Ethylene 
(Flow  Rate  10,000  SCCM,  NO  Cone.  170  ppm.  Water  Absent,  Ethylene  Absent) 
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Figure  14,  Summary  of  NO  Removal  for  Low  Concentration  Data 
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;ure  28.  Ethylene  Concentration  Profiles  in  Wet  and  Dry  Air 
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Figure  30.  Glycolaldehyde  Concentration  Profiles 
(k, =  5E+12) 
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Appendix  A:  Computer  Programs  for  Modeling  Reactions 

This  Program  is  used  to  call  the  ODE  Solver  and  uses  the  conp.inp  file, 
chem.irtp  files  as  input  files. 

This  program  is  for  the  mole  fractions  as  its  input  . . 

PROGRAM  COMP 

Integration  of  adiabatic,  constant  pressure  kinetics  problems 
VERSION  1.2: 

1.  Implement  new  VODE  solver 
*****precision  >  double 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z),  INTEGER { I-N) 

*****END  precision  >  double 
*****preCision  >  single 

IMPLICIT  REAL  (A-H,0“Z),  INTEGER  (I-N) 

*****END  precision  >  single 

PARAMETER  (LENIWK=4000,  LENRWK=4600,  LENCWK=500,  NK=5,  NLMAX^55, 

1  LIN=5,  L0UT=6,  LINKCK-25,  KMAX=50,  ITOL^l,  IOPT=0, 

2  RTOL-1 . OE~6,  ITASK=1,  ATOL=l .OE-15 ) 

DIMENSION  IWORK (LENIWR) ,  RWORK (LENRWK) ,  X(KMAX),  Z (KMAX) 

CHARACTER  CWORK (LENCWK) *  16,  KSYM (KMAX) *16,  LINE*80 
LOGICAL  KERR,  IERR 
EXTERNAL  FUN 

COMMON  /RCONS/  P,  RU 
COMMON  /ICONS/  KK,  NWT,  NH,  NWDOT 

DATA  KERR/. FALSE./,  X/KMAX*0.O/,  KSYM/KMAX* '  V 

WRITE  (LOUT,  15) 

15  FORMAT  ( 

It'  CONP:  CHEMKIN-II  Version  1.2,  Aug.  1992’, 
t****precision  >  double 

C  2/f  DOUBLE  PRECISION') 

C*****END  precision  >  double 
C*****precision  >  single 
C  2/'  SINGLE  PRECISION') 

C*****END  precision  >  single 
C 

C  Initialize  CHEMKIN 

OPEN  (19,  FILE=,'Conp. dat", STATUS=” UNKNOWN ") 

OPEN  {22,  FILE-f,COMP.  DAT1",  STATUS“,,UNKNOWN,, } 

OPEN  (LINKCK,  FORM- ' UNFORMATTED ’ ,  STATUS- 1  UNKNOWN ' , 

1  FILE— ’ chem . bin ' ) 

CALL  CKLEN  (LINKCK,  LOUT,  LENI,  LENR,  LENC) 

CALL  CKINIT  (LENIWK,  LENRWK,  LENCWK,  LINKCK,  LOUT,  IWORK, 

1  RWORK,  CWORK) 

CLOSE  (LINKCK) 

CALL  CKINDX  (IWORK,  RWORK,  MM,  KK,  II,  NFIT) 
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o  o  o  o  o  n  o  o  o 


C 


NEQ 

= 

KK  +  1 

LRW 

22  +  9* NEQ  +  2*NEQ**2 

NVODE 

— 

LENR  +  1 

NWT 

— 

NVODE  +  LRW 

NR 

NWT  +  KK 

NWDOT 

= 

NH  +  KK 

NTOT 

= 

NWDOT +  KK  -  1 

LIW 

= 

30  +  NEQ 

IVODE 

= 

LENI  +  1 

ITOT 

- 

IVODE  +  LIW  -  1 

IF  (KK 

.GT.  KMAX)  THEN 

WRITE  (LOUT,  *) 

1  '  Error... KMAX  too  small... must  be  at  least  ',KK 

KERR  =  .TRUE. 

ENDIF 

C 

IF  (LENRWK  .LT.  NTOT)  THEN 
KERR  =  .TRUE. 

WRITE  (LOUT,  *) 

1  '  Error. . .LENRWK  too  small... must  be  at  least',  NTOT 

ENDIF 
C 

IF  (LENIWK.  .LT.  ITOT)  THEN 
KERR  =  .TRUE. 

WRITE  (LOUT,  *) 

1  ’  Error. . .LENIWK  too  small... must  be  at  least',  ITOT 
ENDIF 
C 

IF  (KERR)  STOP 
C 

CALL  CKSYMS  (CWORK,  LOUT,  KSYM,  IERR} 

IF  (IERR)  KERR  «  .TRUE. 

CALL  CKWT  (I WORK,  RWORK ,  RWORK(NWT)) 

CALL  CKRP  (IWORK,  RWORK,  RU,  RUC,  PATH) 

Pressure  and  temperature 

WRITE  (LOUT,  '(/A)')  ’  ADIABATIC  FIXED  PRESSURE  PROBLEM' 

WRITE  (LOUT,  M/A)')  '  INPUT  PRESSURE  (ATM)  AND  TEMPERATURE (K) ’ 

READ  (LIN,  *)  PA,  T 

WRITE  (LOUT, 7105)  PA,  T 
WRITE  (19,*)  PA,  T 
P  «  PA* PATH 

Initial  non-zero  moles 

40  CONTINUE 
LINE  «  1  ' 

C  WRITE  (LOUT,  M/A)')  ’  INPUT  MOLES  OF  NEXT  SPECIES’ 

READ  (LIN,  '(A)',  END=45)  LINE 
C  WRITE  (LOUT,  MIX,  A)'}  LINE 

ILEN  =  INDEX  (LINE,  M ’ ) 
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c 


IF  {ILEN  .EQ.  1)  GO  TO  40 


I LBN  =  ILEN  -  1 

IF  {ILEN  .LE.  0)  ILEN  =  LEN (LINE) 

IF  { INDEX (LINE { i ILEN) ,  'END')  .EQ.  0)  THEN 
IF  {LINE  ( : ILEN)  .NE.  ’  ' )  THEN 

CALL  CKSNOM  (LINE (s ILEN) ,  1,  LOUT,  KSYM,  KK,  KNUM, 
I  NVAL,  VAL,  I ERR) 

IF  ( I ERR)  THEN 

WRITE  (LOUT,*)  '  Error  reading  moles.,,' 

KERR  -  .TRUE. 

ELSE 

X(KNUM)  -  VAL 

ENDIF 
ENDIF 
GO  TO  40 
ENDIF 

45  CONTINUE 

Final  time  and  print  interval 

WRITE  {LOUT,  '(/A)’)  '  INPUT  FINAL  TIME  AND  DT ' 

READ  {LIN,  *)  TTl,T2,  DT,  III 
WRITE  {LOUT, 7105)  T2 ,  DT,  III 

WRITE  (22,7100)  (KSYM(K),  K=1,KK) 

IF  {KERR)  STOP 

Normalise  the  mole  fractions 

XTOT  »  0.00 
DO  50  K  =  1,  KK 

XTOT  =  XTOT  +  X{K) 

50  CONTINUE 

DO  55  K  =  1,  KK 

X (K)  =  X (K)  /  XTOT 
55  CONTINUE 

WRITE (19,*)  XTOT 

Initial  conditions  and  mass  fractions 

TT1  =0.0 
Z(l)  =  T 

CALL  CKXTY  (X,  I WORK,  RW0RK ,  Z(2)} 

Integration  control  parameters  for  VODE 

TT2  =  TT1 

MF  =  22 
ISTATE=  1 
NLINES=NLMAX  +  1 

Integration  loop 

250  CONTINUE 

IF  (NLINES  .GE.  NLMAX)  THEN 
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c 

C  Print  page  heading 

C 

C  WRITE  (LOUT,  7003) 

C  WRITE  (LOUT,  7100)  (KYSM(K) ,  K=1,KK) 

C  NLINES  *  1 

C 

C  DO  200  K1  =  NK+1,  KK,  NK 

C  WRITE  (LOUT,  7110)  (KSYM.(K)  ( :  10)  ,  K«K1,  MIN (Kl+NK-1,  KK)  ) 

C  NLINES  =  NLINES  +  1 

C  200  CONTINUE 

C  ENDIF 

C 

C  Print  the  solution 
C 

T  =  2(1) 

CALL  CKYTX  (2(2) .  IWORK,  RWORK,  X) 

C 

WRITE  (LOUT,  7105)  TTl,  T,  (X(K),  K=1,KK) 

NLINES  =  NLINES  +  1 

DO  300  K1  »  NK+1,  KK,  NK 

WRITE  (LOUT,  7115)  (X(K),  K=K1,  MIN (Kl+NK-1, KK) ) 

NLINES  =  NLINES  +  1 
300  CONTINUE 

IF  (TT2  . GE.  T2)  THEN 
WRITE (13, 7012 ) (KSYM (K) ,X (K) ,  K=1,KK) 

WRITE (19,*) 'END* 

JJJ=III/2 
QQQ=JJJ*2 

IF(QQQ.EQ.III)  THEN 
XINT=lE-6 
T4-TT2 
TT1=T2 
T2-T2+XINT 
T3=XINT/10 

IF  ( MOD (111,120) . EQ . 0 )  THEN 
WRITE (25, 7105)  TT2 , T, (X (K) ,  K=  1,KK) 

ENDIF 
ELSE 

XINT«16.667E-3 
TTl=T2 
T2=T2+XINT 
T3-XINT/10 
ENDIF 
111=111+1 

WRITE (19,*)  TTl,  T2, T3 ,  III 
STOP 
ELSE 

TT2  »  MIN (TT2  +  DT,  T2) 

ENDIF 


Call  the  differential  equation  solver 
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350  CONTINUE 

C*****precision  >  single 
C  CALL  SVODE 

C*****END  precision  >  single 
C****"* precision  >  double 
CALL  DVODE 

C*****END  precisioii  >  double 

*  (FUN,  NEQ,  Z,  TT1,  TT2 ,  ITOL,  RTOL ,  ATOL,  1TASK, 

1  I STATE,  IOPT,  RWORK(NVGDE)  ,  LRW,  I WORK ( IVODE} , 

2  LIW,  JAC,  MF,  RWORK,  IWORK) 

C 

IF  {I STATE  .LE.  -2)  THEN 
IF  { I STATE  .EQ.  -1)  THEN 
I STATE  «  2 
GO  TO  350 

ELSE 

WRITE  (LOUT,*)  *  ISTATE= 1 , ISTATE 
STOP 

ENDIF 
ENDIF 
GO  TO  250 
C 

C  FORMATS 


7003  FORMAT 
7100  FORMAT 
7105  FORMAT 
7110  FORMAT 
7115  FORMAT 
7012  FORMAT 
END 


( 1 H 1  > 

(3X,  ' T  (SEC)  ' ,  9X,  T TMP (K) 

(36E16.8) 

(26X,  5  (IX, AID)) 

(22X,  10E11.3! 

(  (IX, AlO) , 4X, 1E14 .8) 


9X,  20 (4X, A10) ) 


SUBROUTINE  FUN  (N.  TIME,  Z,  ZP,  RPAR,  IPAR) 
C 


C*****preCision  >  double 

IMPLICIT  DOUBLE  PRECISION (A-H, O-Z ) ,  INTEGER ( I~N) 
C*****END  precision  >  double 
C*****precision  >  single 
C  IMPLICIT  REAL  (A-H, O-Z),  INTEGER (I-N) 

C*****END  precision  >  single 
C 

COMMON  /RCONS/  P,  RU 

COMMON  /ICONS/  KK,  NWT,  NH,  NWDOT 

DIMENSION  Z(*),  ZP<*),  RPAR ( * ) ,  IPAR(*) 

C 

C  Variables  in  Z  are:  2(1)  =  T 
C  Z (K+l )  =  Y (K) 

C 

C  Call  CHEMKIN  subroutines 
C 


CALL  CKRHOY  (P,  Z  (1) ,  Z{2),  IPAR,  RPAR,  RHO) 

CALL  CKCPBS  (Z<1),  2(2},  IPAR,  RPAR,  CPB) 

CALL  CKWYP  (P,  Z(l),  2(2),  IPAR,  RPAR,  RPAR (NWDOT) ) 
CALL  CKHMS  (Z(l),  IPAR,  RPAR,  RPAR (NH) ) 
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Form  governing  equation 


SUM  -  0.0 
DO  100  K  =  1,  KK 

H  =  RPAR {NH  +  K  -  1} 
WDOT  =  RPAR (NWDOT  +  K  -  1) 
WT  =  RPAR (NWT  +  K  -  1) 

ZP (K+l )  -  WDOT  *  WT  /  RHO 
SUM  »  SUM  +  H  *  WDOT  *  WT 
100  CONTINUE 

SP(1)  -  -SUM  /  (RH0*CP8) 

RETURN 

END 
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#  This  program  executes  the  Pulse-on  or  Pulse  off  duration  of  the  Program 
# ! /bin/sh 

#  to  execute:  sh  conp.sh  logname  & 
sh  1>  $  UK  log  2>&1  «  ENDSH 

set  ~x 

led  / s c r / $  LOGNAME 
#mkdir  "${1}$$" 

#cd  ”${1}$$" 

cat  «  EOF  >  makefile 

include  chemmake.h 

OBJS  =  conp.o  cklib.o  vode.o  math.o 
IMPS  =  therm.dat  chem.inp  conp.inp 
OUTS  =  chem.out  chem.bin  conp.out 
EXES  =  chem.exe  conp.exe 

chem.exe:  ckinterp.o 

$ (LINK)  chem.exe  ckinterp.o 

conp.exe:  ${OBJS) 

${LIMK)  conp.exe  ${OBJS) 

EOF 


touch  makefile;  make  chem.exe  conp.exe 
chem.exe;  conp.exe  <  conp.inp  >  conp.out 
cat  fort. 2 5  »  yy2.out  & 
rm  fort. 25 

mv  Conp.dat  conp.inp 
mv  chem.inp  temp . inp 
mv  tchem.inp  chem.inp 
mv  temp . inp  tchem.inp 
leat  conp.out  >>  yy.out  & 
mv  conp.out  temp. out  & 

#f77  “O  average  average. f 
leverage 

cat  tempi. out  >>  temp2.out  & 

ENDSH 


#go  to  user’s  scratch  directory 
Intake  subdirectory  /my run## 

#go  to  /my run## 
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#  This  program  runs  the  whole  simulation  using  different  programs  and 
files 

#f/bin/csh  -f 
set  icount  =  1 
set  iend  *  7201 
while  ($icount  !=  $iend) 
sh  cconp.sh  cconp 
@  ieount++ 
end 

#f77  — o  ave  ave.f 
lave 

f77  ~o  init  init.f 
init 

#cat  INITIAL, DAT  yy2.out  >  yyl.out 
#rm  yy2.out 
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#  This  is  an  example  of  the  input  file. 

#  The  first  line  indicates  the  pressure  in  atm  and  Temperature  in  kelvin. 

#  The  following  lines  indicate  the  initial  mole  fraction  of  different 
species 

#  Finally  the  line  after  the  End  statement  refers  to  initial  time,  final 
time 

#  the  time  interval  and  the  number  of  run  respectively. 


1.00000000000000 

298.000000000000 

02 

0.21 

0 

0. 

n 

0. 

N2 

0.79 

NO 

0.0001 

03 

0. 

E 

0.1E-09 

N02 

0. 

N03 

0. 

N20 

0. 

END 

0  IE-06  1E-Q7 

1 
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Appendix  B:  Senior  Honors  Thesis  -  The  Effects  of  Particles  on  Pulsed  Streamer 

Corona  Discharge  Characteristics 
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Abstract 


Pulsed  streamer  corona  treatment  is  an  advanced  oxidation  technology  using  a 
nonthermal  plasma  that  produces  hydroxyl  radicals,  hydrogen  peroxide,  and  aqueous 
electrons,  all  of  which  react  with  water-borne  organic  contaminants  in  the  removal 
process  (Clements  el.  al.,  1985).  Pulsed  streamer  corona  is  an  electrical  discharge  created 
by  applying  a  continual  series  of  sustainable,  short  lived,  high  voltage  pulses  to  a  small 
diameter  metallic  electrode.  This  discharge  creates  an  electric  field  between  a  point-to- 
plane  electrode  geometry  within  a  batch  reactor  containing  organic  contaminants  in 
deionized  water.  A  rotating  spark  gap  power  supply  generates  this  pulsed  corona,  which 
has  electrical  characteristics  of  short  (200  -  1000  ns)  pulse  width,  fast  (20  -  50  ns)  rise 
time,  and  60  Hz  pulse  frequency 

Recent  experimental  results  have  shown  that  the  addition  of  various  types  of  solid 
particles  to  test  solutions  in  the  aqueous  phase  pulsed  streamer  corona  reactor  affects  the 
physical  nature  of  the  corona  discharge.  The  composition  of  these  solutions  included 
adding  granular  activated  carbon,  (2360  pm  particle  diameter),  powdered  activated 
carbon  (75  -  300  pm  average  particle  diameter),  and  small  glass  beads  (of  sizes:  110- 
180  pm,  53  -  78  pm,  30  -  64  pm,  and  25  and  finer  particle  diameters)  to  deionized  water 
in  a  one  liter  volume  reactor  vessel.  The  discharge  and  ground  electrodes  are  suspended 
within  the  solution.  The  addition  of  the  powdered  carbon  to  the  aqueous  system  showed 
a  change  in  characteristics  of  the  corona  discharge  such  as  increased  numbers  of 
streamers,  an  increase  in  streamer  length,  and  most  notably,  a  large  increase  in  the 
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maximum  applied  voltage  achieved  prior  to  sparkover.  With  the  addition  of  the  glass 
particles,  the  discharge  characteristics  of  the  corona  remained  the  same,  and  the 
sparkover  voltage  increased  with  the  addition  of  the  two  smaller  dimensioned  glass 
spheres. 
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Introduction 


Recently,  there  has  been  an  increased  interest  in  the  efficient  and  cost  effective 
removal  of  organic  contaminants  from  groundwater  and  wastewater.  Hazardous  organic 
contaminants  such  as  phenol,  benzene,  and  PCBs  are  found  in  surface  and  subsoils  and 
also  in  the  treated  sewage  effluent  contained  in  storage  vessels.  Several  sources 
contribute  to  this  problem  including  leaking  petroleum  tanks,  landfill  and  cropland 
runoff,  and  illegal  industrial  waste  dumping.  Federal  and  state  laws  now  mandate  that 
these  contaminants  be  removed  and  that  clean  water  be  returned  to  the  environment. 

Investigations  in  the  field  of  removing  organic  contaminants  from  aqueous 
solutions  have  included  using  advanced  oxidation  technologies  such  as  direct  ozonation, 
electron  beam  bombardment,  UV  photolysis,  ultrasonification,  microwave  plasma 
reduction,  and  gamma  radiation  treatment.  Over  the  past  decade,  another  technology 
involving  the  use  of  a  pulsed  streamer  corona  discharge  in  the  liquid  phase  (Clements,  et 
al.,  1985)  has  been  found  to  be  an  effective  way  to  remove  organics  from  aqueous 
solutions.  This  technique  utilizes  chemical  radicals  produced  from  a  pulsed  streamer 
corona  discharge  that  is  sustained  in  an  aqueous  phase  medium.  Preliminary  studies  at 
the  FAMU/FSU  College  of  Engineering  have  demonstrated  the  complete  degradation  of 
phenol  in  solution  (Sharma,  1993),  and  the  determination  of  the  rates  of  formation  of 
hydrogen  peroxide,  hydroxyl  radicals,  and  aqueous  electrons  (Joshi,  1994). 

Recently,  it  has  been  observed  in  our  laboratory  that  the  addition  of  particles  to  a 
test  solution  affects  the  physical  characteristics  (i.e.,  streamer  length,  intensity,  number  of 
streamers,  and  sparkover  voltage)  of  the  streamer  corona  discharge.  The  principal 
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objective  of  the  present  work  is  to  document  preliminary  experimental  results  of  the 
effects  of  the  addition  of  several  different  types  of  particles  to  the  aqueous  phase  pulsed 
streamer  corona  reactor.  The  particles  used  for  this  investigation  included  1)  75  -  300pm 
diameter  powdered  activated  carbon,  2)  1.40  -  3.35  mm  diameter  granular  activated 
carbon,  3)  1 10  - 180  pm  diameter  glass  beads,  4)  53  -  78  pm  diameter  glass  beads,  5)  30  - 
64  pm  diameter  glass  beads,  and  6)  25  pm  and  finer  diameter  glass  beads.  The 
experiments  conducted  in  this  work  focus  on  determining  how  the  physical  nature  of  the 
corona  discharge  (size  and  number  of  streamers)  is  affected  by  the  addition  of  these 
particles  to  the  pulsed  corona  reactor,  and  how  the  maximum  voltage  before  sparkover  is 
changed  when  the  particles  are  added. 

In  this  report,  a  brief  literature  review  in  chapter  II  will  give  a  background  on  the 
pulsed  corona  process.  In  chapter  III,  a  simplified  model  of  the  process  will  be  given. 

The  experimental  equipment  and  procedure  will  be  discussed  in  chapter  IV.  The  results 
of  a  number  of  experiments  will  be  presented  and  discussed  in  chapter  V.  Chapter  VI 
will  discuss  future  experimental  and  theoretical  work. 
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Literature  Review 


Pulsed  corona  technology  has  been  employed  in  pollution  migration  for  almost  a 
decade  and  a  half.  The  use  of  a  pulsed  streamer  corona  has  proven  to  be  an  effective 
method  of  removing  contaminants  such  as  NOx  and  SOx  from  air  streams  (Masuda, 

1990;  Kamase,  1991).  Clements  et  al.  (1985)  were  the  first  to  adapt  a  pulsed  streamer 
corona  to  an  aqueous  phase  system.  This  group  also  performed  work  on  ozone 
formation,  dye  breakdown,  and  electrical  discharge  characteristics  in  an  aqueous  phase 
pulsed  streamer  corona  system.  The  discharge  consists  of  repetitive  pulses  with  a  short 
pulse  width  (200  ns)  at  an  applied  voltage  of  20  -  45  kV.  The  advantage  of  using  a 
pulsed  power  supply  over  continuous  power  (DC)  is  that  the  energy  input  to  the  system  is 
used  to  produce  reactive  chemical  radicals  while  minimizing  the  heating  of  the  solution. 
The  chemical  radicals  produced  from  the  pulsed  streamer  corona  discharge  in  the 
aqueous  phase  are  given  in  Figure  1 . 

More  recently,  the  aqueous  phase  pulsed  streamer  corona  technology  has  been 
significantly  extended,  and  experimental  and  modeling  aspects  of  this  work  have  been 
documented  in  several  publications  (Sharma  et  al.,  1993;  Joshi  et  al.,  1994).  Many  of  the 
chemical  reactions  involved  in  this  process  have  been  characterized  (Clements  et  al.; 
1985,  Sharma  et  al.,  1993;  Joshi  et  al.,  1994).  These  reactions  include  such  species  as 
hydroxyl  radicals,  hydrogen  peroxide,  and  aqueous  electrons.  The  reaction  and  removal 
rates  of  phenol  have  been  reported  for  different  experimental  conditions  (Sharma  et  al., 
1993).  Experiments  were  carried  out  using  a  point-to-plane  electrode  geometry  within  a 
batch  reactor  filled  with  deionized  water  and  various  additives.  In  all  cases,  at  higher 
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applied  electric  fields  the  rate  of  phenol  removal  increased.  The  addition  of  iron  to 
the  reactor  was  found  to  increase  the  rate  of  phenol  removal  due  to  Fenton's  reaction 
(Joshi,  et  al.,  1994;  Sharma,  et  al.,  1993).  Several  other  additives,  such  as  bubbling 
oxygen  through  the  discharge  electrode,  enhanced  the  rate  of  phenol  removal. 

The  present  work  concerns  the  addition  of  coarse  and  fine  particles  to  the  aqueous 
phase  reactor  to  ascertain  the  particles'  effect  on  the  sparkover  voltage  and  the  nature  of 
the  discharge.  The  addition  of  particles  to  the  pulsed  corona  system  has  not  been 
reported  in  the  literature,  however,  the  addition  of  particles  can  affect  both  the  chemical 
and  physical  aspects  of  the  discharge.  The  preliminary  experiments  presented  in  this 
report  indicate  that  adding  particles  (powdered  activated  carbon,  in  particular)  to  the 
reactor  greatly  affects  the  properties  of  the  corona  discharge.  To  better  understand  this 
phenomenon,  the  corona  discharge  and  the  propagation  of  the  streamers  need  to  be 
analyzed. 

Analysis  of  the  voltage  breakdown  of  a  pulsed  electrode  system  in  aqueous 
solutions  has  been  previously  presented  and  modeled  in  several  different  ways.  One 
method  of  analysis  has  been  through  an  approximation  of  a  shock  wave  (Klimkin  et  al., 
1978).  In  this  method,  shock  waves,  similar  to  ultrasonic  treatments,  are  used  to  estimate 
the  basic  energy  relations  for  the  beginning  stage  of  breakdown  in  the  aqueous  phase. 
Another  mechanism  was  proposed  that,  during  pulse  application,  hydrogen  atoms  become 
excited  in  the  gas  phase,  and  lead  to  the  formation  of  gas  bubbles  (Kuskova,  1983).  This 
formation  has  been  observed  in  previous  work  (Clements  et  al.,  1985).  These  gas 
bubbles  expand  and  have  secondary  reactions  such  as  hydroxyl  radicals  and  peroxides 
(Klimkin,  1990). 
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A  method  proposed  and  developed  by  Neimeyer  et  al.  for  modeling  gas  phase 
dielectric  discharge  has  been  introduced  by  studying  computer  simulations  of  stochastic 
Monte  Carlo  models  (Niemeyer  et  al.,  1984).  To  do  this,  they  solved  the  following 
Laplace  equation: 

V24>=0  (1) 

over  a  circular  lattice  (matrix),  the  center  potential  (discharge  needle  point)  being  zero 
and  the  outer  ring  potential  (ground  plate)  being  one.  Each  time  this  equation  is  solved,  a 
probability  for  the  streamer  propagation  to  an  adjacent  point  in  the  grid  is  given  to  each 
of  the  corresponding  points  on  the  matrix.  Propagation  is  then  achieved  using  a  random 
move  for  each  step.  The  Laplace  equation  is  solved  for  each  iteration.  Based  upon  this 
model,  more  complicated  models  have  been  developed  to  determine  the  effect  of  other 
variables  on  the  nature  of  the  breakdown  (Niemeyer  et  al.;  1989,  Kupershtok,  1992).  The 
addition  of  particles  and  spatial  discontinuities  to  the  grid  developed  by  Neimeyer  et  al. 
(1989)  may  allow  for  a  qualitative  analysis  of  the  effect  of  various  sizes,  shapes,  and 
material  types  of  particles.  Also  note  that  the  material  type  may  also  alter  the  dielectric 
constants. 
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Modeling 

A  simple  model  was  derived  from  a  charge  balance  to  simulate  the  system  used  in 
the  present  work.  Using  Poisson’s  Equation,  with  the  given  boundary  conditions, 


d24>  x  =  0,  4>  =  1 

- =  -p  (2) 

dx 2  x  =  1,  cj>  =  0 


where  4>  is  the  electric  potential,  and  p,  assumed  to  be  a  constant,  is  the  charge  density  of 
the  solution,  we  can  obtain  a  simple  one-dimensional  model  for  the  electrical  discharge  in 
the  system  for  this  work.  The  solution  of  this  equation  is  the  following: 

=  l(p/2)x  +  i](i-x)  0) 

The  graph  of  this  solution  is  shown  in  Figure  2.  In  this  figure,  each  line  represents  a 
constant  charge  density.  From  this  figure,  it  is  noticed  that  as  p  increases,  the  maximum 
electric  potential,  <$>m,  increases.  When  the  graph  of  maximum  potential,  <|>m,  versus 
different  values  of  the  charge  density,  p,  Figure  3  is  obtained.  This  gives  a  relation 
between  the  charge  density  in  the  solution  containing  particles,  and  the  maximum  electric 
potential. 

The  maximum  potential  is  expected  to  simulate  the  voltage  at  which  sparkover 
occurs.  The  charge  density  can  be  related  to  the  surface  area  of  the  particles  in 
suspension  and  to  the  composition  of  the  particles  themselves.  The  following  equations 
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for  field  charging  will  show  how  charge  density  is  related  to  the  surface  area  of  the 

particle  and  the  composition  of  the  particle.  When  a  single  dielectric  particle  of  radius  Rp  ^ 

containing  a  charge  q,  is  placed  in  an  electric  field,  Ew,  in  cylindrical  coordinates,  at  any 
point  outside  the  sphere  is  (Stratton,  1941,  Flagan  and  Seinfeld,  1988): 

q  \  r  Ki-  k2  RP2  1 

V= - +| - + -  |£«,/?pcos0  (4) 

4tt€oK2 r  L  Rp  ki+2ki  r  J 

where  Ki  is  the  dielectric  constant  of  water,  K2  is  the  dielectric  constant  of  the  particle 

added,  and  €o  is  the  permittivity  of  free  space.  The  values  for  K  are  80  for  water,  5.75  for 
carbon,  and  4.42  for  the  glass  beads.  Differentiating  this  equation  in  the  r-component: 

dV  r  Ki  -  K2  r  Rp3 1  q 

Ee= - =  -EooCOSS  I  1  +2 -  |+ -  (5) 

dr  L  K1+2K2/*5  J  4TTr2€oK2 

The  charge  on  the  particle  increases  as  more  charge  migrates  to  the  particle.  Under  the 
conditions  present  in  the  pulsed  corona  reactor,  the  saturation  charge  is  attained  soon 
after  initiation  of  the  discharge.  The  final  charge  on  the  particle  is  represented  by  the 
following: 
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(6) 


r  3ki  i 

q/c  =  I - I  TreoK2£ooZ?P2 

[  Ki+2k2  J 

For  the  limiting  case  of  equation  (6),  where  Ki  >>  K2 ,  such  as  the  system  in  this  project, 
the  following  equation  arises: 

q/c=  3h€oK2£ooZ?p2  (7) 

which  relates  the  field  charge  to  K2,  which  is  the  dielectric  constant  of  the  particle.  In 

equation  (7),  the  given  charge  density  is  linearly  proportional  to  the  surface  area  through 
Dp2,  where  Dp  is  the  particle  diameter.  To  determine  the  charge  density  from  the  field 
charge,  multiply  the  particle  charge  by  the  total  number  of  particles  in  the  solution,  and 
then  divide  this  total  by  the  volume  of  the  solution.  This  gives  a  relationship  between  the 
charge  density,  surface  area  of  the  particle,  and  the  composition  of  the  particle.  These 
charge  densities  can  now  be  applied  to  the  theoretical  models  from  equation  (3). 


Experimental  Equipment  and  Procedure 


The  experimental  apparatus  used  in  this  study  is  shown  in  Figure  4.  The 
apparatus  consists  of  a  550  ml  plexiglass  vessel  used  as  a  reaction  chamber.  This  vessel 
is  closed  at  both  ends  and  is  submerged  in  an  ice  bath  to  approximate  isothermal 
conditions.  The  region  of  electrical  discharge  is  located  at  the  tip  of  a  hollow  hypodermic 
needle,  which  is  connected  to  a  high  voltage  pulsed  power  supply.  The  corona  extends 
from  the  needle  upwards  toward  the  stainless  steel  ground  plate  in  a  point-to-plane 
geometry.  The  separation  distance  between  the  needle  tip  and  ground  plate  is  set  at  5  cm. 
A  magnetic  stirrer  at  the  bottom  of  the  vessel  provides  thorough  mixing  of  the  solution. 
Two  nylon  fittings  on  the  side  of  the  vessel  are  used  to  take  samples  for  analysis. 

The  difference  between  this  system  and  other  corona  discharges  such  as  dc 
discharge,  ac  discharge,  and  long-pulse  (ps  -  ms  pulse  width)  treatments  is  that  the 
electrical  discharge  in  this  reactor  has  a  short  pulse  width  and  a  fast  rise  time.  Pulsed 
power  is  produced  by  a  rotating  spark  gap  high  voltage  pulsed  power  supply,  capable  of 
supplying  the  reactor  discharge  electrode  with  high  voltage  (20  -  45  kV),  short  duration, 
(200  -  1000  ns),  fast  rise  time  (20-50  ns),  repetitive  (60  Hz ),  electrical  pulses.  The 
pulsed  power  supply  used  for  this  project  is  shown  in  Figure  5.  The  source  input  is  a 
high  voltage  transformer  supplying  a  rectified  AC  power  to  the  pulse  forming  circuit. 

This  high  voltage  goes  through  a  series  of  current-limiting  resistors  and  a  set  of  diodes 
that  remove  the  negative  part  of  the  sine  wave.  The  half-wave  rectified  voltage  then 
charges  an  array  of  ceramic  "doorknob"  capacitors  which  discharge  into  the  reactor  load 
whenever  the  rotating  spark  gap  closes. 
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For  each  set  of  experimental  conditions,  three  or  more  trials  were  made  to  obtain 
appropriate  statistical  interpretations  of  the  data.  The  initial  experimental  trials  for  this 
work  consisted  of  control  runs  of  550  ml  of  deionized  water  only  in  the  pulsed  corona 
reactor.  Two  different  controls  were  run,  during  which  measurements  of  the  voltage  and 
current  between  the  high  voltage  and  ground  electrodes  were  made  and  qualitative 
observations  of  the  length  and  quantity  of  the  streamers  were  made.  The  first  control 
utilized  a  fully  exposed  discharge  electrode  needle,  and  the  second  employed  a  needle 
covered  with  a  plastic  covering  along  the  needle  shaft,  leaving  only  the  tip  exposed  in 
order  to  simulate  a  more  true  point-to-plane  geometry. 

The  first  set  of  experiments  testing  the  effects  of  adding  particles  to  the  reactor 
solution  were  performed  using  powdered  activated  carbon  (50  to  200  mesh,  300  -  75  pm 
particle  size)  at  concentrations  of  0.5  g,  1.0  g,  2.0  g,  3.0  g,  and  5.0  g  per  550  ml  of 
deionized  water.  The  carbon  was  added  to  the  stirred  water  in  the  vessel  before  the  trials 
began  to  obtain  a  well-mixed  solution  when  the  experiment  started.  The  second  set  of 
experiments  used  a  larger  type  of  particle,  granular  activated  carbon  (6  to  14  mesh,  3.3  - 
1 .4  mm  particle  size).  These  experiments  were  run  with  concentrations  of  1 .0  g,  2.0  g, 
and  5.0  g  of  carbon  per  550  ml  of  deionized  water.  Finally,  the  last  set  of  experiments 
were  performed  with  four  different  sizes  of  spherical  silica  glass  beads,  1 10  -  180  pm 
diameter,  53  -  78  pm  diameter,  30  -  60  pm  diameter,  and  25  pm  and  finer  diameter  glass 
beads.  Again,  each  of  the  sizes  were  evaluated  at  concentrations  of  1.0  g,  and  2.0  g,(and 
5.0  g  at  the  1 10  - 180  pm  diameter  and  53  -  78  pm  diameter  sizes  only)  per  550  ml 
deionized  water. 
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Once  these  test  solutions  were  made,  samples  were  taken  to  measure  the  initial 
temperature,  pH,  and  conductivity.  To  begin  a  trial,  a  test  solution  was  introduced  into 
the  reactor,  and  the  pulsed  power  supply  was  engaged  at  an  initial  applied  voltage  of  10 
kV.  The  voltage  was  then  increased  at  increments  of  2.5  kV  until  sparkover  was  reached. 
Sparkover  occurs  when  the  applied  voltage  is  high  enough  to  cause  a  spark  to  bridge  the 
two  electrodes.  At  each  voltage  setting,  the  average  peak  current  was  read  using  an 
ammeter,  and  visual  observations  were  made  of  the  size,  quality,  and  quantity  of  the 
streamers.  Subsequent  to  sparkover,  samples  were  again  taken  to  test  the  temperature, 
pH,  and  conductivity  of  the  solution. 

Also,  aqueous  solutions  of  potassium  iodide  (KI)  were  made  at  concentrations  of 
0.0301  M.  These  solutions  contained  no  particles  but  had  high  ionic  content.  These 
experiments  were  done  to  compare  the  affects  of  a  solution  with  a  higher  conductivity 
due  to  a  dissolved  salt,  to  that  of  solid  particles  in  suspension.  Each  of  these  solutions 
was  tested  before  and  after  treatment  for  temperature,  pH,  and  conductivity.  The 
experiments  with  the  KI  started  with  an  applied  voltage  of  10  kV  and  increased 
incrementally  by  2.5  kV  until  approximately  80  kV. 
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Experimental  Results  and  Discussion 


Current  vs.  voltage  plots  for  the  two  initial  experiments  with  no  particles  in  the 
system  are  shown  in  Figure  6.  The  streamers  observed  emanating  from  the  discharge 
electrode  appeared  higher  in  number  for  the  control  run  with  the  covered  needle  electrode 
when  compared  to  the  run  with  the  uncovered  needle.  Sparkover  was  reached  at  42.5  kV 
for  the  uncovered  needle  control  and  sparkover  was  at  40  kV  for  the  covered  needle 
control.  These  observations  are  consistent  with  those  found  in  previous  work  (Sharma,  et 
al.,  1993).  All  further  experiments  were  conducted  with  the  covered  needle,  and  these 
experiments  were  compared  to  the  data  and  observations  of  the  control  run  with  the 
needle  shaft  covered  in  plastic.  Once  again,  the  covered  needle  gives  the  system  a  truer 
point-to-plane  geometry  than  using  the  uncovered  needle. 

The  first  set  of  experiments  testing  particles  in  solution  were  performed  by  adding 
powdered  activated  carbon  at  concentrations  of  0.5  g,  1.0  g,  2.0  g,  3.0  g,  and  5.0  g  to  550 
ml  of  deionized  water  in  the  pulsed  corona  reactor.  The  current  vs.  applied  voltage  data 
are  presented  in  Figure  7.  In  all  five  trials  with  the  particles  present,  the  current 
measurement  at  the  corona  onset  voltage  of  20  kV  was  higher  than  that  for  the  no-particle 
control  trial.  All  of  the  trials  with  particles  produced  a  greater  number  of  streamers  and 
longer  streamers  at  a  given  voltage  when  compared  to  the  control  with  no  particles. 

These  determinations  are  qualitative  by  observation  and  no  exact  measurements  were 
made.  As  the  concentration  of  the  powdered  carbon  increased  from  0.5  g  to  5.0  g„  the 
quantity  and  length  of  streamers  also  increased  for  equal  voltage  applied. 

The  most  significant  observations  were  the  voltage  differences  at  which  sparkover 
occurred.  Sparkover  is  very  distinguishable  in  that  there  is  a  very  bright  streamer 
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connecting  the  two  electrodes  accompanied  by  a  loud  noise.  Table  1  shows  the  voltage  at 
sparkover  for  the  six  different  powdered  carbon  particle  concentrations. 


Table  1.  Carbon  Particle  Concentration  With  Sparkover  Voltages 


Concentration  of  Carbon  Particles 

per  550  ml  DI  Wager  (g) 

0 

0.5  g 

1.0  g 

2.0  g 

3.0  g 

5.0  g 

Sparkover  Voltage  (kV) 

40  kV 

45  kV 

82.5  kV 

>100  kV 

This  table  shows  that  as  the  concentration  of  the  carbon  particles  is  increased,  the 
applied  voltage  at  which  the  system  achieves  sparkover  also  increases.  This  phenomena 
is  significant  because  the  higher  the  electric  field  (i.e.  applied  voltage)  the  more  chemical 
radicals  produced  and  therefore  a  higher  contaminant  removal.  Treatment  at  or  above 
sparkover  essentially  short  circuits  the  system  therefore  not  producing  as  many  of  the 
chemical  radicals  that  the  pulsed  streamer  corona  creates. 

These  results  show  a  large  increase  in  the  maximum  applied  voltage  achieved 
before  sparkover  as  the  concentration  of  the  powdered  activated  carbon  is  raised.  The 
results  also  give  a  graph  (Figure  7)  that  shows  a  linear  increase  in  the  applied  voltage  at 
sparkover  as  the  amount  of  carbon  increases.  In  previous  work  in  this  laboratory  using 
only  deionized  water  in  the  reactor,  the  maximum  voltage  prior  to  sparkover  was  40  kV 
for  the  same  electrode  geometry  without  particles.  With  the  addition  of  the  carbon 
particles,  higher  voltages  can  be  achieved,  which  may  lead  to  higher  conversion  rates  and 
faster  degradation  of  organic  contaminants  in  the  reactor.  The  solution  with  the  5  g  of 
carbon  never  reached  sparkover  because  the  upper  voltage  limit  of  the  power  supply  is 
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100  kV.  On  Figure  7,  the  graphs  show  a  clear  picture  of  two  different  slopes  for  each  of 
the  lines  representing  current,  a  small  slope  initially,  and  then  at  a  given  point,  a  larger 
slope  occurs. 

The  experimental  results  for  the  sparkover  voltage  vs.  the  total  particle  surface 
area  are  shown  in  Figure  8.  A  best  fit  of  these  points  gives  a  straight  line.  This  linear 
trend  corresponds  with  that  predicted  from  the  theoretical  background  for  large  p  shown 
in  Figure  3.  For  these  quantities  of  carbon  particles,  the  total  surface  area  of  the  particle 
samples  are  very  high  and  therefore  the  range  of  the  charge  density  corresponds  with  the 
straight  portion  of  the  slope  in  Figure  3.  There  is  a  need  to  relate  the  concentrations  of 
the  carbon  particles  to  the  charge  density. 

The  second  type  of  particle  tested  in  the  corona  reactor  was  the  granular  activated 
carbon.  These  experiments  used  samples  of  550  ml  deionized  water  containing 
concentrations  of  1 .0  g,  2.0  g,  and  5.0  g  of  granular  activated  carbon.  Figure  9  shows  the 
relationship  between  the  current  and  voltage  for  these  four  granular  carbon  concentration 
experiments.  There  were  no  significant  differences  between  the  granular  activated 
carbon  runs  and  the  control  except  for  the  5  g  concentration  experiments.  At  5  g  of 
carbon,  the  measured  current  was  slightly  higher  and  the  sparkover  increased  to  42.5  kV 
from  40  kV  at  all  other  concentrations.  There  were  no  observed  differences  in  streamer 
quantity  or  length  between  any  of  the  granular  carbon  runs.  A  possible  reason  for  the 
lack  of  change  in  the  properties  of  the  corona  might  be  that  the  granular  activated  carbon 
particles  were  difficult  to  suspend  in  solution,  and  vigorous  stirring  caused  most  of  the 
particles  to  be  accelerated  outward  to  the  reactor  margin,  away  from  the  middle  of  the 
solution  where  the  high  voltage  discharge  needle  is  located. 
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The  next  set  of  experiments  were  done  with  two  different  sizes  of  glass  beads, 

1 10  - 180  pm  diameter  and  53  -  78  pm  diameter.  Again,  each  of  the  sizes  was  evaluated 
at  concentrations  of  1.0  g,  2.0  g,  and  5.0  g  per  550  ml  deionized  water  in  the  pulsed 
corona  reactor.  The  effects  of  the  particles  on  the  current  vs.  voltage  characteristics  are 
shown  in  Figures  10  and  1 1 .  For  fixed  voltages,  there  was  an  increase  in  current  for 
experiments  with  the  smaller  sized  glass  beads  and  less  of  a  change  for  the  experiments 
with  larger  sized  glass  beads  when  compared  to  the  no-particle  control.  The  sparkover 
voltage  was  37.5  kV  for  all  six  of  the  concentrations  and  sizes  of  glass  beads.  The  most 
notable  observation  for  all  of  these  trials  was  that  the  streamer  size  decreased,  but  the 
number  of  streamers  increased  and  they  seemed  to  be  much  more  intense  and  brighter. 

At  these  concentrations,  the  current  varied  for  each  of  the  runs.  This  is  why  the  statistical 
analysis  is  not  shown  on  any  of  the  graphs.  A  reason  for  this  problem  is  that  the  ammeter 
might  have  been  damaged  from  overuse  and  this  could  be  the  reason  why  the  results  were 
sporadic. 

Also,  the  smaller  glass  particles  that  were  looked  at  were  the  30  -  64  pm  diameter 
glass  beads  and  the  25  pm  diameter  and  finer  glass  beads.  The  results  for  the  current  vs. 
voltage  graphs  is  shown  in  Figures  12  -  15.  These  figures  show  the  same  trend  as  the 
previous  glass  beads  in  that  the  current  increases  as  the  voltage  increases.  The  25  pm  and 
finer  diameter  glass  beads  showed  the  closest  resemblance  to  the  powdered  activated 
carbon  in  that  the  sparkover  voltage  increased  with  increasing  concentrations  of  the  glass 
beads.  When  the  glass  beads  of  this  diameter  are  suspended  by  the  stirring  bar,  the 
solution  becomes  very  cloudy  therefore  visual  observations  of  the  corona  discharge  were 
not  made. 
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Sparkover  voltage  plotted  against  the  surface  area  for  all  of  the  glass  particles  is 
given  in  Figure  16.  The  glass  particles  had  very  little  total  surface  area  when  compared 
to  the  carbon  particles  due  to  the  carbon  particles  being  very  porous  and  the  glass  is  non- 
porous.  Due  to  this  small  surface  area,  this  graph  may  correspond  to  the  section  of  the 
graph  of  charge  density  vs.  surface  area  in  the  theoretical  background  (Figure  3).  The 
graph  in  Figure  3  first  shows  a  flat  section  and  then  an  increase  in  the  maximum  electrical 
potential  at  small  values  of  the  charge  density. 

An  analysis  of  the  charge  density  is  shown  in  Figure  17.  This  is  shown  by 
plotting  the  sparkover  voltage  vs.  the  charge  density  of  the  various  concentrations  of 
glass  beads  in  solution.  Charge  density  was  found  by  first  finding  the  charge  density  on 
each  particle  using  the  field  charging  theory  in  equation  (6),  multiplying  this  number  by 
the  total  number  of  particles,  and  then  dividing  by  the  volume  of  the  reactor.  The  result 
of  this  analysis  is  a  slope  that  is  similar  to  that  of  the  predicted  theoretical  graph  of  the 
maximum  electric  potential  versus  the  charge  density  as  shown  in  Figure  3.  The  graphs 
look  similar  over  the  initial  regions  where  the  slope  is  flat  and  then  starts  to  increase 
corresponding  to  small  values  in  p.  The  solutions  for  the  carbon  were  not  analyzed  in 
this  manner  because  the  field  charging  equation  is  modeled  for  a  system  for  solid 
spherical  particles  in  an  electric  field.  The  carbon  particles  are  not  spherical  and  were 
very  porous.  The  actual  total  surface  area  is  106  times  greater  than  the  outer  surface  area. 
An  equation  needs  to  be  devised  that  will  account  for  this  great  surface  area. 

Finally,  the  results  of  the  trials  with  the  potassium  iodide  solutions,  the  current  vs. 
voltage  plots  are  given  in  Figure  18.  These  experiments  were  done  to  compare  the  affects 
of  solution  with  a  higher  conductivity  due  to  a  dissolved  salt  to  that  of  solid  particles  in 
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suspension.  The  graph  in  Figure  18  shows  an  increase  in  current  with  an  increase  in 
applied  voltage.  This  system  did  not  achieve  sparkover  even  when  the  applied  voltage 
reached  80  kV.  The  conductivity  of  this  solution  was  104  times  greater  than  the  solutions 
with  the  carbon  and  the  glass  particles.  The  streamers  in  this  solution  were  much  smaller 
than  DI  water  alone,  but  the  streamers  were  much  more  intense. 

The  current  to  voltage  relationships  for  all  of  the  particles  had  different  slopes. 
The  KI  solutions  had  the  smallest  slopes,  the  glass  beads  had  a  slightly  higher  slope,  DI 
water  was  greater  than  the  glass  beads,  and  the  powdered  activated  carbon  particles  had 
the  highest  slope.  These  carbon  particles  were  also  the  only  solution  that  had  two 
different  slopes.  The  slope  of  the  current  for  these  carbon  particles  started  at  one  slope 
and  then  suddenly  increased  the  slope  of  the  current.  The  cause  of  these  slopes  is 
unknown  at  this  time.  The  conductivity  of  the  solution  might  have  a  function  in  the 
different  slopes  of  the  solutions. 

From  these  results,  the  particle  type  having  the  greatest  effect  on  the  sparkover 
voltage  in  the  reactor  was  the  powdered  activated  carbon.  Addition  of  powdered  carbon 
to  the  aqueous  solution  produced  more  streamers,  longer  streamers,  and  increased 
sparkover  voltage.  With  these  characteristics,  it  is  possible  that  the  removal  of  organics 
such  as  phenols  would  be  accelerated  greatly  since  the  production  of  hydroxyl  radicals, 
aqueous  electrons,  and  hydrogen  peroxide  increases  with  increases  in  the  applied 
discharge  voltage.  This  would  lead  to  decreased  exposure  time  and  therefore  increased 
efficiency. 
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Future  Work 


These  results  bring  up  many  new  questions  on  the  addition  of  particles  to 
deionized  water  in  the  pulsed  streamer  corona  reaction  vessel.  More  investigation  of  the 
physical  aspects  of  the  process  as  well  as  the  chemical  aspects  of  the  process  should  be 
investigated. 

Further  development  of  the  model  presented  in  this  work  is  needed.  The  least 
understood  part  of  this  model  is  the  definition  of  the  maximum  electric  potential.  It  is 
predicted  that  the  maximum  potential  is  related  to  the  point  at  which  sparkover  is 
achieved,  but  the  relation  is  not  well  understood.  Also,  other  kinds  of  particles  such  as 
copper  or  photocatalysts  should  be  looked  at  to  aid  in  this  prediction.  The  copper  will 
give  a  non-porous  conducting  particle  and  the  photocatalysts  could  promote  surface 
reactions. 

To  this  point,  chemical  aspects  of  the  effect  of  adding  particles,  such  as  the  rates 
of  formation  of  the  reactive  species,  and  the  rate  of  phenol  breakdown,  has  not  been 
investigated.  Pulsed  corona  treatments  of  contaminated  water  containing  particles  should 
be  fully  analyzed.  Carbon  is  known  to  adsorb  small  aromatic  organics  and  carbon  has 
shown  to  promote  surface  reactions  (Vidic,  Suldan,  and  Brenner,  1993.  Abuzald  and 
Nakhla,  1994.)  The  kinetics  of  the  formation  of  the  different  reactive  species  should  be 
reanalyzed  for  the  addition  of  the  carbon  and  glass  particles.  Also,  the  adsorption  of  the 
contaminants  and  subsequent  surface  reactions  should  also  be  noted. 
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Applied  Voltage  (kV) 


Current  vs  Applied  Voltage  for  Powdered  Activated  Carbon 
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Current  vs  Applied  Voltage  for  110  - 180  pm  Glass  Beads 
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Current  vs  Applied  Voltage  for  53  -  78  pm  Diameter  Glass  Beads 
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Experimental  Results  for  Sparkover  Voltage  vs  Charge  Density 
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